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1 Introduction 



Over the past two decades Cosmology has increasingly become a precision science. That the 
Universe is expanding was an astonishing discovery. Now we know its details to unprecedented 
precision. An expanding Universe also implied an extremely compact state in the past, and 
therefore very high temperature. The Particle Physics forces which can now be explored only in 
accelerator laboratories were in free play in the remote past. Thus the observation of the oldest 
remnants in the Universe amounts to looking at the results of a Particle Physics experiment under 
natural conditions. 

In these notes we present a selection of topics, each section approximately amounting to one 
lecture. We begin with a brief recapitulation of General Relativity, and the Standard Model of 
Cosmology. The study of Cosmology requires General Relativity to be applied only under a highly 
symmetric situation and therefore it is possible to recast the essentials as Three Laws of Cosmology. 
The study of very early Universe brings us squarely into the domain of Quantized Field Theory 
at given temperature. Intermediate metastable phases through which the Universe passed require 
an understanding of the effective potential of the field theory in a thermal equilibrium. This 
formalism is developed in some detail. 

The remainder of the notes discuss important signatures of the remote past. These include : (i) 
inflation, (ii) density perturbations leading to galaxy formation, (iii) study of hot and cold relics 
decoupled from the remaining constituents, some of which can be candidates for Dark Matter, (iv) 
finally the baryon asymmetry of the Universe. As we shall see each of these has a strong bearing 
on Particle Physics and is being subjected to ever more precise observations. 

1.1 General Theory of Relativity : A Recap 

Special Theory of Relativity captures the kinematics of space-time observations. On the other 
hand, General Theory of Relativity is a dynamical theory, which extends the Newtonian law of 
gravity to make it consistent with Special Relativity. In this sense it is not a "generalization" 
of Relativity but rather, a theory of Gravity on par with Maxwell's theory of Electromagnetism. 
It is nevertheless a very special kind of theory because of the Principle of Equivalence. The 
equivalence of gravitational and inertial masses ensures that in a given gravitational field, all test 
particles would follow a trajectory decided only by their initial velocities, regardless of their mass. 
This makes it possible to think of the resulting trajectory as a property of the spacetime itself. 
This was the motivation for introducing methods of Differential Geometry, the mathematics of 
curved spaces for the description of Gravity. Due to this "grand unification" of space and time 
into a dynamically varying set we shall use the convention of writing the word spacetime without 
hyphenation as adopted in rcf. [1]. 

Throughout these notes we use the convention h = c = 1 and the sign convention ds 2 = 
dt 2 — \dx\ for the spacetime interval in Special Relativity. The principle of General Covariance 
states that a given gravitational field is described by a metric tensor g^ v a function of spacetime 
variables collectively written as x^ = x, t. Gravity modifies the spacetime interval to the general 
quadratic form g^dx^dx" ', where the summation convention on same indices is assumed. The 
trajectories of test particles are simply the shortest possible paths in this spacetime, determined 
by the metric tensor through the geodesic equation 



These symbols are not tensors but determine the covariant derivative much the same way that 
the electromagnetic potentials which are themselves not gauge invariant determine the minimal 
coupling of charged particles to electromagnetic fields. 



d 2 x^ 




where the Christoffel symbols T% are given by 
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The equations which determine the gravitational field, i.e., the tensor itself are the Einstein 
Equations, 

G>„ — A<7 M „ = R^ v — -^9fj,vR — A<7 M „ — SirGT^v 

where T^ v is the energy momentum tensor and the Ricci tensor R^ u and the scalar curvature R 
are the contracted forms of the fourth rank tensor the Riemann curvature, given by 

^uap — U <x L v0 U P L va T 1 cot v(3 1 afj 1 va 

Rpv = R^xv 
R = gTR^ 

The tensor G^ v is called the Einstein tensor and has the elegant property that its covariant 
derivative vanishes. The last term on the left hand side is called the Cosmological term since its 
effects are not seen on any small scale, even galactic scales. It can be consistently introduced 
into the equations provided A is a constant. Since the covariant derivative Dpg^ = 0, the 
covariant derivative of this term also vanishes. This fact is matched on the right hand side of 
Einstein's equation by the vanishing of the covariant derivative of the energy-momentum tensor. 
The vanishing of the Einstein tensor follows from the Bianchi identities in Differential Geometry. 
The geometric significance of the identities is that given a spacetime domain, they express the 
statement "the boundary of a boundary is zero". We leave it to the reader to pursue ref [1] to 
understand its details. Thus a geometric principle implies the covariant conservation of energy- 
momentum tensor, a physical quantity. But it has to be noted that covariant conservation does 
not imply a conserved charge the way it happens in flat spacetime with divergence of a four- vector. 
But if there are 3-dimensional regions on whose 2-dimensional boundaries gravity is very weak, it 
does imply conservation of total mass-energy in the given volume. 

There are quite a few subtleties concerning the implications of General Relativity and the 
conditions under which it supercedes Newtonian gravity. We present here a few "True or False" 
statements for the reader to think over and discuss with peers or teachers. Starting points to 
answers are given in Appendix, sec. 1121 

1.1.1 True or False 

1. Curved spacetime is a necessity of GR due to the underlying Special Relativity principle. 

2. The invariance of the equations of physics under arbitrary reparameterisation of spacetime 
is the essential new content of GR. 

3. The notion of energy density becomes meaningless in GR 

4. The notion of total energy becomes meaningless in GR 

5. Points where some of the metric coefficients vanish are unphysical 

6. Points where some of the metric coefficients diverge are unphysical 

7. Points where any components of curvature tensor diverge are unphysical 

8. Newtonian gravity is insufficient to describe an expanding Universe and GR is required. 

1.2 The Standard Model of Cosmology 

Here we summarise the broadest features of our current understanding of the Universe based on a 
vast variety of data and modeling. We summarise size, age and contents of the Universe as follows 
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• There is a strong indication that the Universe is homogeneous and isotropic if we probe it 
at sufficiently large scales, such as the scale of clusters of galaxies. The typical scale size 
is 10 Megaparsec (Mpc) and larger. At the scale of several tens of Mpc the distribution of 
galaxies is homogeneous. 

It is at present not known whether the Universe is finite in size or infinite. A finite Universe 
would be curved as a three dimensional manifold. An infinite universe could also show 
curvature. At present we do not see any signs of such curvature. 

• Secondly we believe that the Universe has been expanding monotonically for a finite time in 
the past. This gives a finite age of about 13.7 billion years to our Universe. 

What was " before" this time is not possible to understand within the framework of classical 
General Relativity. But Newton's gravitational constant suggests a fundamental mass scale, 
called the Planck scale Mpi — G -1 / 2 sa 1.2 x 10 19 GeV and corresponding extremely small 
time scale, 10 -44 sec. We expect that Quantum theory of Gravity should take over at that 
scale. Unfortunately that theory has not yet been worked out due to insufficient data. 

• Finally it has been possible to map the current contents of the Universe to a reasonable 
accuracy. More about them below. The contents of the Universe can be divided into three 
types, 

1. Radiation and matter as known in Particle Physics 

2. Dark Matter, one or more species of particles, at least one of which is necessarily non- 
relativistic and contributing significantly to the total energy density of the Universe 
today. These particles have no Standard Model interactions. 

3. Dark Energy, the largest contributor to the energy density balance of the present Uni- 
verse, a form of energy which does not seem to fit any known conventional fields or 
particles. It could be the discovery of a non-zero cosmological constant. But if more 
precise observations show its contribution to be changing with time, it can be modelled 
as a relativistic continuum which possesses negative pressure. 

Of the contents of type 1, there are approximately 400 photons per cc and 10 -7 protons 
per cc on the average. Compared to Avogadro number available on earth, this is a very sparse 
Universe. Of these the major contributor to energy density is baryonic matter. This constitutes 
stars, galaxies and large Hydrogen clouds out of which more galaxies and stars are continuously 
in formation. The other component is the Cosmic Microwave Background Radiation (CMBR), 
the gas of photons left over after neutral Hydrogen first formed. Its contribution to total energy 
density is relatively insignificant. But its precision study by experiments such as the Wilkinson 
Microwave Anisotropy Probe ( WMAP) is providing us with a very detailed information of how 
these photons have been created and what are the ingredients they have encountered on their way 
from their origin to our probes. 

There probably are other exotic forms of energy-matter not yet discovered. Two principle 
candidates are topological defects such as cosmic strings, and the axion. We shall not be able to 
discuss these here. Axions are almost a part of Standard Model though very weakly interacting, 
and could potentially be Dark Matter candidates. 

Aside from these current parameters of the Universe there is a reason to believe that the 
Universe passed through one or more critical phases before arriving at the vast size and great 
age it presently has. This critical phase in its development is called Inflation. It is expected to 
have occurred in remote past at extremely high energies, perhaps in the Planck era itself. What 
is interesting is that the fluctuations in energy density which finally became galaxies could have 
originated as quantum effects during that era. Thus we would be staring at the results of quantum 
physics in the very early Universe whenever we see galaxies in the sky. 

A quantitative summary of present observables of the Universe is given at the end of sec. [3] 
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1.3 The Standard Model of Particle Physics 

We assume that the reader is familiar with the Standard Model of Particle Physics. Appendix B of 
Kolb and Turner contains a review. Cosmology has a dual role to play in our understanding of the 
fundamental forces. It is presenting us with the need to make further extensions of the Standard 
Model and is also providing evidence to complete our picture of Particle Physics. For example 
whether Dark Matter emerges from an extension of the Standard Model is a challenge to model 
building. On the other hand axions expected from QCD may perhaps get verified in Astroparticle 
physics experiments and may have played a significant role in the history of the Cosmos. 

Study of Cosmology also sharpens some of the long recognized problems and provides fresh 
perspectives and fresh challenges. Symmetry breaking by Higgs mechanism in Standard Model 
(and its extensions the Grand Unified models) causes hierarchy problem. But it also implies a 
cosmological constant far larger than observed. We hope that the two problems have a common 
solution. Despite the conceptual problems with the QFT of scalar fields, inflation is best modeled 
by a scalar field. Similarly, consider Dark Energy which is almost like a cosmological constant of 
a much smaller value than Particle Physics scales. This has finally been confirmed over the past 
decade. Again many models seem to rely on unusual dynamics of a scalar field to explain this 
phenomena. We hope that supersymmetry or superstring theory will provide natural candidates 
for such scalar fields without the attendant QFT problems. 

2 Friedmann-Robertson- Walker Metrics 

Cosmology began to emerge as a science after the construction of reflection telescopes of 100 
to 200 inch diameter in the USA at the turn of 1900. When Doppler shifts of Hydrogen lines 
of about twenty nearby galaxies could be measured it was observed that they were almost all 
redshifts. Edwin Hubble proposed a linear law relating redshift and distance. Then the data 
could be understood as a universal expansion. Over the last 75 years this fact has been further 
sharpened, with more than 10 million galaxies observed and cataloged. 

It is reasonable to believe that we are not in a particularly special galaxy. So it is reasonable to 
assume that the expansion is uniform, i.e. observer on any other galaxy would also see all galaxies 
receding from him or her at the same rate. This reasoning allows us to construct a simple solution 
to Einstein's equations. We assume that the Universe corresponds to a solution in which aside 
from the overall expansion as a function of time, the spacetime is homogeneous and isotropic. This 
would of course be true only for the class of observers who are drifting along with the galaxies in a 
systematic expansion. ( An intergalactic spaceship moving at a high speed would see the galactic 
distribution quite differently). We characterize this coordinate system as one in which (1) there 
is a preferred time coordinate such that at a given instant of this time, (2) the distribution of 
galaxies is homogeneous and isotropic. 

It should be noted that this is a statement about the symmetry of the sought after solution 
of the Einstein equations. The symmetries restrict the boundary conditions under which the 
equations are to be solved and in this way condition the answer. In older literature the existence 
of such symmetries of the observed Universe was referred to as Cosmological Principle, asserted in 
the absence of any substantial data. Over the years with accumulation of data we realize that it is 
not a principle of physics, but a useful model suggested by observations, similar to the assumption 
that the earth is a spheroid. 

The assumptions of homogeneity and isotropy made above are very strong and we have a very 
simple class of metric tensors left as possible solutions. They can be characterized by spacetime 
interval of the form 



1 This time coordinate is not unique. We can define a new time coordinate t(t) as an arbitrary smooth function 
of the old time coordinate t. What is unique is that a time axis is singled out from the spacetime continuum. What 
units, (including time dependent ones), we use to measure the time is arbitrary. 
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The only dynamical degree of freedom left is the scale factor R(t). Further, there are three 
possibilities distinguished by whether the spacelike hypersurface at a given time is flat and infinite 
( Newtonian idea), or compact like a ball of dimension 3 with a given curvature, (generalization 
of the 2 dimensional shell of a sphere in usual 3 dimensional Euclidean space), or unbounded and 
with constant negative curvature, a possibility more difficult to visualize. These three possibilities 
correspond to the parameter k = 0, or k = +1, or k = —1. The cases k = — 1 and k = 1 also have 
representations which make their geometry more explicit 

( dt 2 - R 2 (t ) { d X 2 + sin 2 x (d0 2 + sin 2 6d(f> 2 ) } k = 1 
ds 2 = <^ 

[ ^-^(^l^+sin^x^+sin 2 ^ 2 )} k = -1 

The time coordinate t we have used above is a particular choice and is called comoving time. 
An alternative time coordinate rj is given by 

dt 

dv = W) 

ds 2 = R 2 ( r] ){drf-^—-dr 2 -r 2 ^ 2 6d(j> 2 
I 1 — kr z 



Its advantage is that for k = it makes the metric conformally equivalent to flat (Minkowski ) 
space. 

2.1 Cosmological redshift 

In this and the next subsection we identify the precise definitions of redshift and cosmological 
distances to understand Hubble Law in its general form. 

The observed redshift of light is similar to Doppler shift, but we would like to think that it 
arises because spacetime itself is changing and not due to relative motion in a conventional sense. 
In other words, in cosmological context the redshift should be understood as arising from the 
fact that time elapses at a different rate at the epoch of emission t\ from that at the epoch to 
of observation. Since light follows a geodesic obeying ds 2 = 0, it is possible for us to define the 
quantity f{r\) which is a dimensionless measure of the separation between emission point r = 
and the observation point r = n 

dr f to dt 

o (l-kr 2 ) 1/2 ~ J tl R(t) 



/(n) S / 



where the second equality uses ds 2 = 0. Now the same f(r\) is valid for a light signal emitted a 
little later at t = t\ + St\ and received at a corresponding later time to + St - 



Equivalently, 



to+St a dt 

t!+6t! R(t) 



rti+st ldt ^ rto+st dt 

iti R J to 



R(t) 



Sti Sto 
R(h) = R(t ) 
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or 



8tx _ Ai 
5t A 



R(t ) 



It is convenient to define the redshift z, originally so denned because it would always be small, as 
given by 

Aq R{to) 



1 + -a7 



2.2 Luminosity Distance 

Defining a measure of spacelike distances is tricky in Cosmology because physical separations 
between comoving objects are not static and therefore lack an operational meaning. Since distances 
are measured effectively by observing light received, we define luminosity distance dh by 



dl = 



L — > absolute luminosity 
AttF — > observed flux 



If the metric were frozen to its value at time t this would have been the same as in flat space, 
R(t ) 2 r 2 with n the coordinate distance travelled by light. Due to expansion effects, we need 
additional factors of 1 + z, once for reduction in energy due to redshift and once due to delay in 
the observation of the signal 

d 2 L =R 2 (t )r 2 (l + z) 2 

We now introduce measures H for the first derivative, representing the Hubble parameter and 
a dimesionless measure qo for the second derivative, traditionally deceleration, by expanding the 
scale factor as 

■^L = l + H (t-t a )- l -q Hl{t-t f + .... 



R(t ) _ R(t ) 

.: z = H (t -t) + (l + f) Hi (to - tf + 
(to-t) = H^ (*-(l + f) Z 2 -{ 



-R 



RHq 



so that 



We now use the quantity f(r) introduced in the discussion of redshift of light, which for the three 
different geometries works out to be 



fin) 




\r\ 3 



Substitutc (t — ti) 



R(t ) 



(to - h) + ^H (t - hf + ... 



n 



l 



R (to) Ho 



H d L = z+^(l-qo)z 2 + ... 
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The last relation expresses the connection between observed luminosity of distance <1l of galaxies 
and their redshift z, incorporating the curvature effects arising from the expansion of the Universe. 
Extensive data on cLl and z gathered from galaxy surveys such as 2-degree Field Galactic Redshift 
Survey ( 2dF GRS) and 6dF GRS can be fitted with this equation to determine cosmological 
parameters H and q . 

3 The Three Laws of Cosmology 

It is possible to discuss the cosmological solution without recourse to the full Einstein equations. 
After all a comprehensive framework like electromagnctism was discovered only as a synthesis of 
a variety of laws applicable to specific situations. Specifically, Coulomb's law is most useful for 
static pointlike charges. Similarly given that the Universe has a highly symmetric distribution of 
matter, allows us to express the physics as three laws applicable to Cosmology, these being 

1. Evolution of the scale factor : The evolution of R(t) introduced above is governed by 

{ ilY k A 8tt 
[r) + &-3=T Gp 

2. Generalized thermodynamic relation : The energy-matter source of the gravitational field 
obeys generalization of the first law of thermodynamics dU = —pdV, 

|<^ )+ 4 (#)-„ 

To put this to use we need the equation of state p = p(p). For most purposes it boils down 
to a relation 

p = wp 

with w a constant. 

3. Entropy is conserved : (except at critical thresholds) 

d , m d ( R 3 . \ „ 
dt^=dt{-T^ + ^)=° 

The law 2 can be used to solve for p as a function of R with w a given parameter. Then this 
p can be substituted in law 1 to solve for R(t) the cosmological scale factor. 

As for Law 3, in the cosmological context we speak of entropy density s. Thus the above law 
applies to the combination sR 3 . Further, non-relativistic matter does not contribute significantly 
to entropy while for radiation s oc T 3 . Hence we get the rule of thumb 

S = sR 3 oc T(t) 3 R(t) 3 = constant 

Thus Law 3 provides the relation of the temperature and the scale factor. However, critical 
threshold events are expected to have occured in the early Universe, due to processes going out of 
equilibrium or due to phase transitions. The constant on the right hand side of R(t)T{t) equation 
has to be then re-set by calculating the entropy produced in such events. 

This formulation is sufficient for studying most of Cosmology. However if we are familiar with 
Einstein's theory the above laws can be derived systematically. We need to calculate the Einstein 
tensor. The components of the Ricci tensor and the scalar curvature R in terms of the scale factor 
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R(t) are 



Roa — - 



4 



„ R R 2 k \ 

^ = - 6 [r + & + & 9ij 



The resulting Einstein equations contain only two non-trivial equations, the Goo component and 
the Gii component where i = 1, 2, 3 is any of the space components. 

n R fit\ k n „ 

n : 2 r + [r) +w = - BnGp 

It turns out that the ii equation can be obtained by differentiating the 00 equation and using 
the thermodynamic relation 

d (pR 3 ) = -pd (R 3 ) 
Hence we only state the second law rather than the ii equation. 

3.1 Example : Friedman Universe 

An early model due to A. A. Friedman (1922) considers a k — 1 universe with pressureless dust, 
i.e., p = and with A = 0. Then according to Law 2, 

Now let t\ be a reference time so that for any other time t, p(t)R 3 (t) = p(ti)R 3 (ti) Then according 
to Law 1, 

J_fdR\ 2 k_ 8tt R\ 



it 2 ^" 



R 
where 

Rmax = ~^-GpiR 1 

This equation is first order in time, but non-linear. It can be solved by expressing both t and R 
in terms of another parameter r\. The solution is 

R {V) = ^Rmax{l - COST?) t(f]) = ^R m ax(v - sm7 ?) 

It results in a shape called cycloid. The first order equation thus solved can also be thought of as 
a particle in potential V(x) = — x "^ ax with total energy — 1. 

Exercise : ( Friedman 1924) Solve for the scale factor of a A = universe with pressureless 
dust but with negative constant curvature, i.e., k = — 1. 
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Age t 

Hubble constant Hq 
Parameter h 
Tcmb 

n 

O h 2 

"all matter"- 

n b h 2 
n A 

w (=p/p) 



13.7±0.2G. yr, 
h x lOOkm/s/Mpc 
0.71 or 0.73 
2.725 x 10VK 
1.02 ± 0.02 
0.120to0.135 
0.022 

0.72 ± 0.04, 
-0.97 ±0.08 



3.2 Parameters of the Universe 

We now summarise the observable parameters of the Universe as available from several different 
data sets. But first we introduce some standard conventions. First we rewrite the first law 
including the cosmological constant, 



k A 8ttG 
H + i^-3 = — P 

Suppose have a way of measuring each of the individual terms in the above equation. Then with 
all values plugged in, the left hand side must balance the right hand side. Our knowledge of 
the Hubble constant Hq is considerably more accurate than our knowledge of the average energy 
density of the Universe. We express the various contributions in the above equation as fractions 
of the contribution of the Hubble term. This can of course be done at any epoch. Thus, dividing 
out by H 2 , 

1 + J^2 = "A + fl P 

where we define the fractions 

A _8nGp 

In detail, due to several different identifiable contributions to p from baryons, photons and Dark 
Matter, wc identify individual constributions again as fractions of the corresponding Hubble term 
as fif,, fl 7 , and Qdm respectively. Either the sum of the various fi's must add up to unity or the 
k =^ 0. In tabic 13.21 we list the current values of the parameters. 

Determination of Hubble constant has several problems of calibration. It is customary to treat 
its value as an undetermined factor h times a convenient value 100 km/s/Mpc which sets the scale 
of the expansion rate. It is customary to state many of the parameters with factors of h included 
since they are determined conditionally. At present h 2 « 1/2. The Dark Energy component seems 
to behave almost like a cosmological constant and hence its contribution is given the subscript A. 
n being close to unity signifies that we seem to have measured all contributions to energy density 
needed to account for the observed Hubble constant. The parameter w exactly —1 corresponds to 
cosmological constant. Current data are best fitted with this value of w and assumption of cold, 
i.e., non-rclativistic Dark Matter. This Friedman model is referred to as the A-CDM model. 



4 The Big Bang Universe 

By its nature as a purely attractive force, gravity does not generically allow static solutions. Since 
the Universe is isotropic and homogeneous now it is reasonable to assume that that is the way 
it has been for as far back as we can extrapolate. Such an extrapolation however will require 
the Universe to have passed through phases of extremely high densities and pressures, phases 
where various microscopic forces become unscreened long range forces, again subject to overall 
homogeneity and isotropy of the medium. This is the essential picture of a Big Bang Universe, 
extrapolated as far back as the Quantum Gravity era about which we can not say much at the 
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present state of knowledge. However the intervening epochs encode into the medium important 
imprints of the possible High Energy forces that may have been operative. Thus the early Universe 
is an interesting laboratory for Elementary Particle Physics. 



4.1 Thermodynamic Relations 

For a massless single relativistic species, 



30 

e(3) 



gT for bosons 



gT 3 for bosons 



x \ for fermions 



X — for fermions 
4 



where g is the spin degeneracy factor. In the early Universe, a particle species with mass much 
less than the ambient temperature behaves just like a relativistic particlcl. More generally, we 
introduce an effective degeneracy g* for a relativistic gas which is a mixture of various species, 



relat 



4 o^relat 



30 5 * T =3P 



9* = 



El Tj 
9. 



T 



Bose 
species 



7 



(I 

3 

Fermi 
species 



As T continues to drop species with rest masses 3> T become nonrelativistic and stop con- 
tributing to the above. For such species, 



non— relat. 



mini 



m t (rcj(ti)fl 3 (ti)) 
R 3 (t) 



non— relat 



Must" 



4.2 Isentropic Expansion 

We assume particle physics time scales to remain far shorter than expansion time scale -ff(i) -1 . 
Thus we have approximate equilibrium throughout and we have been using equilibrium thermody- 
namics. Specifically entropy is conserved. While the usual term for such a system is adiabatic, the 
term in the sense of being isolated from "other" systems does not apply to the Universe and we 
shall refer directly to fact that entropy is conserved. With no other system with which to exchange 
energy, it can be shown for the primordial contents of the Universe, that (Kolb & Turner, Section 
3.4) 



T 



= dS 



2 Massless vector bosons have one degree of freedom less than massive vector bosons. So this has to be accounted 
for in some way. In the Standard Model of Particle Physics, masses arise from spontaneous symmetry breaking 
which disappears at sufficiently high temperature and the longitudinal vector boson degree of freedom is recovered 
as additional scalar modes of the Higgs at high temperature. 
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we define s = £ = ^t 2 which is dominated by contribution of relativistic particles 



V T 

2TT 2 

s = —g* s T 



with 



TA 3 7v^ /T 7 N ' 



, / j 
Bose Fermi 

Note that this is a different definition of the effective number of degrees of freedom than in the 
case of energy density p. 

4.3 Temperature Thresholds 

While the conditions of approximate equilibrium and isentropic expansion hold for the most part, 
crucial energy thresholds in microphysics alter the conditions in the early Universe from time to 
time and leave behind imprints of these important epochs. Based on the scales of energy involved, 
we present here a short list of important epochs in the history of the early Universe 
T Relevant species g* 

T < MeV 3z/s; photon ( 7) 3.36 (see reason below) 

1 MeV to 100 MeV 3z/s; 7; e+e" 10.75 

> 300 GeV 8 gluons, 4 electroweak gauge 106.75 

bosons; quarks & leptons (3 gen- 
erations), Higgs doublet 
Above 300 GeV scale, no known particles are non-relativistic, and we have the relations 

1/2 T 2 
H = 1.66.9 ' 



M P , 

t = 0.30/ gi> 2 ^~ I 



: >-\ln ( T 2 
" MeV 



sec 



where 

G = 1 



{Mpif (1.2211 x 10 19 GeU) 2 
4.4 Photon Decoupling and Recombination 

We now see in greater detail how one can learn about such thresholds, with the example of 
photons. As the primordial medium cools, at some epoch, neutral hydrogen forms and photons 
undergo only elastic (Thomson) scattering from that point onwards. Finally photons decouple 
from matter when their mean free path grows larger than Hubble distance if -1 . 

T 7 = n e <TT < H 

Here <jt = 6.65 x 10 _25 cto 2 is the Thomson cross-section. (Verify that in eV units, ctt = 1-6 x 
10~ 4 (McV)~ 2 ). Note that we should distinguish this event from the process of "recombination", 
which is the formation of neutral Hydorgen. In this case the competition is between cross section 
for ionisation and the expansion rate of the Universe. By contrast the above relation determines 
the epoch of decoupling of photons, also sometimes called "the surface of last scattering" . Thus the 
process of recombination ends with residual ionisation which is the n e required above to determine 
decoupling. 
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4.4.1 Saha Equation for Ionization Fraction 

To put above condition to use, we need n e as a function of time or temperature and then H at the 
same epoch. n e should be determined from detailed treatment of non-equilibrium processes using 
Boltzmann equations, which we shall take up later. However, utilizing various conservation laws, 
we can obtain a relationship between the physical quantities of interest, as was done by Saha first 
in the context of solar corona. 

We introduce the number densities nn, n pi n e of neutral Hydrogen, protons and neutrons 
respectively. The amount of Helium formed is relatively small, njy e — 0.1n p and is ignored. 
Charge neutrality requires n e — n p , and Baryon number conservation requires ns = nn + n p . 
From approximate thermodynamic equilibrium, we expect these densities to be determined by the 
Boltzmann law, where i stands for H, p or e : 

with c/i degeneracy factors, m,i the relevant masses and ^ the relevant chemical potentials respec- 
tively. Due to chemical equilibrium, 

m = He + Mp 

Thus we can obtain a relation 

9H fm e T\- 3/2 /m H \ 3/2 ( B\ 
n H = n p n e — — exp — 

9p9e V V 2n J \ m pJ \ T J 

where B denotes the Hydrogen binding energy, B — m p + m e — m# = 13.6 eV. Finally we focus 
on the fraction X e = ^ of charged baryons relative to total baryon number, 




Figure 1 : Ionization fraction as a function of photon temperature expressed in eV 
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Fig. \T\ shows variation of A| 9 as function of the photon temperature T, using the value of 
ns/n-f = 6 x 1CP 10 as known from Big Bang Nucleosynthesis calculations and WMAP data. We 
see that the ionization fraction reduces to 5% when T w 0.29eV~ 350K. These are our estimates 
of the ionisation fraction (set arbitrarily at 5%) and the temperature of the last scattering. 

Using this estimate we can now calculate the Hubble parameter in the inequality which de- 
termines decoupling. We assume that the Universe was already matter dominated by this epoch. 
Then the scale factor obeys 

R(t) = ^ 

where to and Rq are a particular epoch and the corresponding scale factor, chosen here to be the 
current values signified by subscript 0. Using this we try to determine the H at decoupling epoch 
in terms of Hq. 



R M Ho \PoJ \n B o 
Combine this with the relation 

Pbo m p n B Q 

J l r _ " 



Pcrito 3H§/8irG 

where m p is the proton mass. Then the condition for photon decoupling becomes 

To \ 3 ^ 2 87r m v 



T J 3 £IbHq 

Substituting ut expressed in eV units computed above, and X^ q of 5% and the corresponding 
temperature found from the graph of Saha formula, we can check that the temperature of the 
photons decoupled at that epoch should today be 

T ~ (32X) x (ft s ) 2/3 

Although very crude, this estimate is pretty close to the order of magnitude of the temperature of 
the residual photons today. We have thus traced a possible way Alpher and Gamow could antici- 
pate the presence of residual radiation from the Big Bang at approximately such a temeperature. 



5 Phase Transitions with Quantum Fields in the Early Uni- 
verse 

We have been considering an expanding Universe, but its expansion rate is so slow that for most 
of its history, it is quasi-static as far as Particle Physics processes are concerned. Under this 
assumption, we can think of a thermal equilibrium for substantial periods on the microscopic 
scale, and build in the slow evolution of the scale factor as a significant but separate effect. 

Quantized fields which are systems of large numbers of degrees of freedom display different 
collective behavior, with qualitative changes in ground state as a function of temperature. The 
technique for studying the ground state behavior at non-zero temperature is a modification of the 
process of calculating Green functions using path integral method in functional formalism. 



5.1 Legendre Transform 

As a warm up for the method to be used, consider a magnetized spin system with spin degrees 
of freedom s(x) as a function of position, and Ti. the Hamiltonian of the system. If we want to 
determine the equilibrium value of magnetization, we first introduce an external field H , in the 
presence of which, the Helmholtz Free Energy F is given by 
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where [3 = 1/T, (same as l/(kT) in our units) and Vs denotes functional integration over s(x). 
Now the quantity of interest is 

/1 r) f)F 
(IlWt)) = __ logZ = -_ 

We now introduce a function of M itself, the Gibbs free energy, through the Legendre transform 

G(M) = F + MH 

so that by inverse transform, 

dM 

Now H being an auxiliary field has to be set to zero. Thus the last equation can be now read as 
follows. The equilibrium value of M can be found by minimizing the function G(M) with respect 
to M. In other words, for studying the collective properties of the ground state, G is the more 
suitable object than H. 

5.2 Effective action and effective potential 

In Quantum Field Theory, a similar formalism can be set up to study the collective behavior of a 
bosonic field </>. It is possible in analogy with the above, to define a functional T of an argument 
suggestively called 4> c i designating the c-number or the classical value of the field. Analgous to 
the external field H above, an auxiliary external current J(x) is introduced. Then 



Z[J] = e- tE ^ = J £V>exp i J d A x (£[<f>] + j(f>) 



Then we obtain the relations 



5E _ .SlogZ 
Jl _ 1 SJ(x) 



-(Sl\<Kx)\Sl). 

4>cl{x) 



Therefore, let 



so that 



r[ct> cl } = -E[j}- J d 4 yJ(y)My) 

^-T [U = -J{x) 

0<Pcl 

Thus the quantum system can now be studied by focusing on a classical field 4> c i , whose dynamics is 
determined by minimizing the functional T. The auxiliary current J is set zero at this stage. This 
is exactly as in classical mechanics, minimizing the action for finding Euler-Lagrange equations. 
The functional T is therefore called the effective action. If it can be calculated, it captures the 
entire quantum dynamics of the field 4> expressed in terms of the classical function <f> c i. 

Calculating T can be an impossible project. A standard simplification is to demand a highly 
symmetric solution. If we are looking for the properties of a physical system which is homogeneous 
and in its ground state, we need the collective behavior of ^ in a state which is both space and time 
translation invariant. In this case 4> c i{x) — > 4> c i, d^4> c i = and T becomes an ordinary function 
(rather than a functional) of (f) c i. It is now advisable to factor out the spacetime volume to define 

Vejf (<f>ci) = -(y^r^ci] 

so that the ground state is given by one of the solutions of 

d 



-V eff (M=0 
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In general V e ff can have several extrema. The minimum with lowest value of V e ff, if it is unique, 
characterises the ground state while the other minima are possible metastable states of the system. 
A more interesting case arises when the lowest energy minimum is not unique. In a quantum system 
with finite number of degrees of freedom, this would not result in any ambiguity. The possibility 
of tunneling between the supposed equivalent vacua determines a unique ground state - as in the 
example of ammonia molecule. But in an infinite dimensional system, such as a field system, we 
shall see that such tunneling becomes prohibitive, and the energetically equivalent vacua can exist 
simultaneously as possible ground states. 

When we speak of several local minima of V e ff and therefore maxima separating them we 
are faced with a point of principle. Any function defined as a Legendre Transform can be shown 
to be intrinsically convex, i.e., it can have no maxima, only minima. The maxima suggested 
by above extremization process have to be replaced by a construction invented by Maxwell for 
thermodynamic equilibria. Consider two minima 4>\ and 4>2 separated by a maximum 03 as shown 
in [21 We ignore the part of the graph containing point 03 as of no physical significance, and 
introduce a parameter x which permits continuous interpolation between the two minima. Over 
the domain intermediate between 0i and 02, we introduce a parameter x and redefine 4>ci to be 

4>d = X(pi + (1 - x)4>2 < x < 1 

The assumption is that the actual state of the system is no longer translation invariant and is an 




Figure 2: The effective potential of a system with local minimum at 0i and 02- The mostly 
concave segment with the point 03 is treated as unphysical and replaced by the dashed line. 

admixture of phases with either a value 0i or a value 02 . The above redefinition is interpreted 
as a value of averaged over such regions. Then if the physical value of C ; is characterised by 
parameter value x, we define the corresponding value of V e ff 

V e a f J (M = xV ef/ (0l) + (1 - x)V eff (4> 2 ) 

These formulae will not be of direct relevance to us. However we should remember that while the 
minima 0i and 02 represent possible physical situations, the extremum does not. We shall see 
that the admixture phase arises when the system begins to tunnel from one vacuum to another 
due to energetic considerations. We shall develop the formalism for the tunneling process and will 
be more interested in the tunneling rate which in turn is indicative of how quickly the transition 
from a false vacuum (a local minimum) to a true vacuum (absolute minimum) can be completed. 
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5.3 Computing V e ff 

A very simple example of a non-trivial effective potential is given by a complex scalar field theory 
with a non-trivial self-coupling 

Note that the mass-squared is negative and so the kinetic terms without the <j> 4 interaction would 
imply tachyonic excitations. The correct field theoretic interpretation is to work around a stable 
minimum such as <f> = fi. Thus we define <p(x) = fi + 4>{x) and quantize only the <j>(x) degrees of 
freedom. By substituting the redefined <j> into the above Lagrangian we find that the field <fi has 
a mass-squared +n 2 and can be treated perturbatively. 

In this example the polynomial (j> — fi 2 ) 2 serves as the zeroth order effective potential. The 
values \ie %a for real values of alpha between and 2tt are all permissible vacuum expectation values 
(VEVs) as determined from minimizing this polynomial. These values are modified in Quantum 
Field Theory, however in this simple example their effect will only be to shift the value [i of \<$>\ 
by corrections of order H. This example is the well known Goldstone mechanism of symmetry 
breaking. Here the symmetry breaking is indicated by the classical V itself. 

There are important exceptions to this, where quantum corrections become important and they 
must be computed. One is where the above Lagrangian is modified as follows. Remove the mass 
term altogether, but couple the massless field to a £7(1) gauge field, in other words massless scalar 
QED. The classical minimum of the potential (f>) 2 is <j> — 0. It was shown by S. Coleman and 
E. Weinberg that after the one-loop quantum effects are taken into account, the minimum of the 
scalar potential shifts away from zero. The importance of this result is that the corrected vacuum 
does not respect U{1) gauge invariance. The symmetry breakdown is not encoded in classical V 
and has to be deduced from the V e f / . 

Another important example is effects of a thermal equilibrium. Consider the complex scalar of 
our first example, with mass-squared negative at classical level. Suppose we couple this field to a 
gauge field. This would lead to spontaneous but classically determinable breakdown of the gauge 
invariance according to Higgs-Kibblc mechanism. If we now include temperature corrections, then 
in the high temperature limit the minimum of the scalar field Lagrangian in fact shifts to <fi = 0, 
thus restoring the symmetry! This is the case important to the early Universe where Universe 
cools slowly from very high values of temperature, possibly close to the Planck scale. Since Higgs 
mechanism is the chief device for constructing Grand Unified Theories (GUTs) we find that gauge 
symmetries are effectively restored at high temperatures and non-trivial vacua come into play only 
at lower temperatures. 

We will be unable to discuss the formalism for computing the effective potential in any de- 
tail. However here we shall briefly recapitulate the recipe which reduces the problem to that of 
computing Feynman diagrams. Recall that an important consequence of quantum corrections can 
be change in the vacuum expectation value of a scalar field. We shall assume that the vacuum 
continues to be translation invariant, so that this shifted value is a constant. Anticipating this, in 
the Lagrangian we shift the field value 4>(x) = <p c i + r/(x) where the significance of the subscript 
cl becomes clear due to the reasoning given above. However note that at his stage 4> c i is only a 
parameter and the V e ff is computed as a function of it. Now in the presence of an external source 
J(x), 

r d A x{C + J4>) = fd A x{C + Jct>d)+ fd 4 xr)(^ + J\ 



Here onwards we assume that the </> c ; is chosen to satisfy 

77 =~J(x) 
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Note that this <p c i is still dependent on the external function J(x) and hence adjustable. Then in 
the path integral formula for the generating functional, we can carry out a saddle point evaluation 
of the gaussiarjf) integral around the extremum defined by the choice of <f> c i just made. 



exp 



Urjexp Jc {<j) cl ) + J<j) c i + \ j V^" vj 
i £(4>d) + J4>d 



x ( det 



-1/2 



where the Z% denotes all the terms of order higher in the variations of C with respect to 
perturbative interpretation hese higher derivatives 

8 n C 



In 



S 2 C 



are treated as vertices, while 



is used as propagator. The main correction to one-loop order however can be determined directly 
from the determinant resulting from saddle point integration. While calculating iS-matrix elements 
this determinant is only an overall constant and of no significance. But in the effective potential 
formalism it provides the main correction. It also requires proof to know that the determinant 
provides the leading correction and that everything inside Z2 represents higher order quantum 
corrections. This will not be pursued here. We now quote an example. 



5.3.1 An example 

Consider the theory of two real scalar fields (f>i and (f>2 

. i 

with i = 1,2 and /1 2 > 0. The latter condition means that the classical minimum of the theory 
is not at fa = but at any of the values defined by \{(t>\ + f^) = /J 2 - A possible minimum is 
at (</>i,02) = (^/\/A, 0). If we shift the fields by choosing (fa, (^2) = {4>d + Vi( x )i V^{ x )) we S e t 
propagators for the two real scalar fields r]i with mass-squares given by m\ = 3A0 2 ; — ^ 2 and 

Evaluation of the determinants of the inverse propagators requires a series of mathematical 
tricks. 

log det (a 2 + rn 2 ) = Tr log (d 2 + m 2 ) 

= y^log (-fc 2 + to 2 ) 

k 

The determinant actually has a diagrammatic interpretation as was explained by Coleman and 
Weinberg. Consider a single closed loop formed by joining together n massless propagators, joined 
together by consecutive mass insertions. Here we are treating treating mass as an interaction 
for convenience. This loop has in front of it from perturbation theory rules due to n mass 

3 Actually "pseudo"-gaussian due to the presence of i, but equivalently after choosing to work with euclidian 
path integral. 



i i 
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insertions. But there are (n — 1)! ways of making these identical mass isertions. This makes the 
contribution of this loop weighed by — . Further because the propagators are identical, nth term 
is nth power of first term. Thus summing all the terms with single loop but all possible mass 
insertions amounts to a log series. The next important argument is that indeed the perturbative 
expansion is an expansion ordered by the number of loops. By including all contributions at one 
loop, we have captured the leading quantum correction. 

We now turn to evaluation of this using dimensional regularization prescription 



/d d k 
-t log (-fc 2 + to 2 ) = 
(2ir) d V ; 



(27T) 



Therefore, 



Veff (<Pd) 



1 2/2 i ^ j.4 
<t>cl + J<t>cl 



(4tt)<V2 ( TO 2)-«V2 



i r (-rf/2) 

"2 (47r) d / 2 



(A</> 2 ; -u 2 ) d/2 + (3A^- M 2 ) 



d/2 



where the <5 M and S\ represent finite parts. Now 

r(2-d/2) 1 (2 



(Air) d / 2 (to 2 ) 



2\2-d/2 



(4tt) 2 



— y + log 4-7T — log TO' 



MS scheme 



(47T) 



log 



TO 
M 2 



where a reference mass scale M has to be introduced. Therefore, with superscript (1) signifying 
one-loop correction, 



V 



(i) 



eff 



1 1 

4 (4tt) 2 



(A^_ u 2 ) 2 (log(A^-n 2 ) /M 2 -3/2) 



+ (3A^ - u 2 ) 2 (log [(3A0 2 , - u 2 ) /M 2 ] - 3/2) 

A rule of thumb summary of this example is that the one-loop correction to the effective potential 
is (1/647T 2 ) m^jjln(m e ff /M) 2 where m 2 ^ contains (j) 2 cl and M is a reference mass scale required 
by renormalization. Further, it can be shown that the modification to the tree level (classical) 
vacuum expectation value comes only from the to 4 ln(m./M) 2 term of the field direction in which 
the vacuum is already shifted, <pi in the present example. 



5.4 Temperature corrections to 1-loop 

In the early universe setting we are faced with doing Field Theory in a thermal bath, also referred 
to as "finite temperature field theory". With some clever tricks and exploiting the analogy of 
the field theory generating functional with the partition function for a thermal ensemble, one 
can reduce this problem also to that of calculating an effective potential. The key modification 
introduced is to combine the time integration of the generating functional and the multiplicative 
—0 (inverse temperature) occurring in the partition function into an imaginary time integral 
— Jq dr. Further, the trace involved in the thermal averaging can be shown to be equivalent to 
periodicity in imaginary time of period (3. For bosonic fields one is lead to periodic boundary 
condition and for fcrmionic fields one is lead to anti-periodic boundary condition. Thus the usual 
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propagator is replaced by an imaginary time propagator of appropriate periodicity. For a scalar 
field of mass m, thermal propagator A T is given by 



k°=2mn/0 ' 



-e 



-ik-(x-y) _ 



I 



We now quote the result for the temperature dependent effective potential V T . We focus only 
on the <pi degree of freedom of the previous example and drop the subscript 1, 



K//[«?y = V ef f [<f> c i] + / dxx 2 \n 







1 — exp — I x 



,2 



771 
y2 



1/2N 



with m (</> c /) = — A* + 3A0 c; 

and 14/ / to one-loop order is as obtained in the previous subsection. In the high temperature limit 
T 3> (pel we can determine the leading effects of temperature by expanding the above expression 
to find 

V?ff = V ef f + ±T><&-^T* + ... 

From thermodynamic point of view this 14// represents the Gibbs free energy. Entropy density 
pressure and the usual energy density are given by 

dV T f( 

s[4>d] = j^ 1 ; whiles -V T [<^] 

Ot 

andp[0 d ] = V^f+Ts^d] 

The resulting graphs of 14// are plotted in fig. [31 In plotting these, the term T 4 which is the 
usual thermodynamic contribution, but which is independent of the flicld <p c i is subtracted. 

We can summarise the main results of this subsection retaining a single scalar degree of freedom 
(j) for which, 



with minima at <j± = ±y/ fJ, 2 /X 

The leading effect of the temperature correction is to add a term T 2 <j) 2 cl so that 

K T //=f-k + k 2 W7A<^ + 



2^ 8 r c A 



This is extremized at the values 



Further, note that 



//i 2 -AT 2 /4 
±\l- : — and 6 d =0 



d 2 V 2 A 

U c =o— -A* + T 1 



As a result we see that the curvature can change sign at the trivial extremum <f> c i — at a critical 
temperature T c = 2/_t/\/A. 
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Figure 3: Temperature dependent effective potential plotted to show its dependence on <j) c i for 
various values of temperature T. T c denotes the temperature below which the trivial minimum is 
unstable. 

6 First order Phase transition 

At the end of the last section we saw that the minimum at (f> = can turn into a maximum as 
the temperature drops below the critical value T c . This effect is felt simultaneously throughout 
the system and a smooth transition to the newly available minimum with <f> ^ ensues. The 
expectation value of <fi is called the order parameter of the phase transition, and in this case where 
it changes smoothly from one value to another over the entire medium is called a Second Order 
phase transition. 

But there can also be cases where there are several minima of the free energy function but 
separated by an energy barrier. Thus the system can end up being in a phase which is a local 
minimum, called the false vacuum, with another phase of lower free energy, called the true vacuum, 
available but not yet accessed. If the barrier is not too high, thermal fluctuations can cause the 
system to relax to the phase of lower free energy. The probability for the system to make the 
transition is expressed per unit volume per unit time and has the typical form 

T = ^exp{-B} 

The expressions A and B are dependent upon the system under consideration. The presence of B 
reminds us of the Boltzmann type suppression that should occur if the system has to overcome an 
energy barrier in the process of making the transition. In Quantum Field Theory there are also 
quantum mechanical fluctuations which assist this process. We observe that the formula above is 
also of the type of WKB transition rate in Quantum Mechanics. Indeed, we have partly thermal 
fluctuations and partly tunneling effects responsible for this kind of transition. A convenient for- 
malism exists for estimating the combined effects using the thermal effective potential introduced 
in the previous section. 

A transition of this type does not occur simultaneously over the entire medium. It is charac- 
terized by spontaneous occurrance of small regions which tunnel or fluctuate to the true vacuum. 
Such regions of spontaneously nucleated true vacuum are called "bubbles" and are enclosed from 
the false vacuum by a thin boundary called the "wall" . Since the enclosed phase is energetically 
favorable, such bubbles begin to expand, as soon as they are formed, into the false vacuum. Over 
time such bubbles keep expanding, with additional bubbles continuing to nucleate, and as the 
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bubbles meet, they merge, eventually completing the transition of the entire medium to the true 
vacuum. Such a transition where the order parameter phi has to change abruptly from one value 
to another for the transition to proceed is called a First Order phase transition (FOPT). 

6.1 Tunneling 

At first we shall consider tunneling for a field system only at T — 0. An elegant formalism has been 
developed which gives the probability per unit volume per unit time for the formation of bubbles 
of true vacuum of a given size. Consider a system depicted in fig. H which has a local minimum at 
value cf>i, chosen to be the origin for convenience. There are other configurations of same energy, 
such as (f>2 separated by a barrier and not themselves local minima. If this is an ordinary quantum 
mechanical system of one variable and the initial value of is 0i , the system is unstable towards 
tunneling to the point 02 and subsequently evolving according to usual dynamics. If this were 

v A 




Figure 4: A system which has a local minimum at 0i and which is unstable towards tunneling to 
a point 4>2 of equal energy 

point particle mechanics, the formula for the transition amplitude from the state \(j>\) to the state 
|02) is given in the Heisenberg picture and in the path integral formulation as 




where T is a time interval and the action S on the right hand side has the same range of time 
integration. Instead of evaluating this directly, we make two observations. Firstly, if we asked for 
the amplitude for the state |0i) to evolve into itself after time T, it would involve contributions 
also from paths that access the state \4>2}- Thus if we inserted a complete set of states on the 
left hand side above at an intermediate time, say, time T/2 (we justify the T/2 later), among the 
many contributions there would also occur the term 

(0i|e-^ T / 2 |0 2 )(0 2 |e-^ T / 2 |0 1 ) 

Correspondingly on the right hand side we would have contribution from paths that start at the 
point 0i and end at 0i, but after reaching 02 somewhere along the trajectory. 

The second point is more interesting. Actually the presence of |02), a state of equal energy, 
makes |0i) unstable. Hence the contribution such as (0i \e~i HT ^ 2 102) actually makes the total 
amplitude for 0i returning to 0i smaller than unity in magnitude. This happens only if the 
evolution operator e~^ HT ^ 2 somehow departs from being of unit magnitude, i.e., its exponent 
becomes real negative rather than pure imaginary. 
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Thus if we look, not for the entire amplitude, but only for the part where the exponent becomes 
effectively imaginary then that part of the sum over intermediate states actually indirectly gives 
the transition amplitude {4>2\e~^- HT / 2 \<f)i), the one we started out to look for (aside from factor 
1/2 in time). In the limit T becomes infinite, all the contributions with real negative exponents 
will go to zero. The leading contribution is the term with the smallest exponent. On the right 
hand side this means that among all the paths that start from one vacuum, sample the other and 
return, the one that minimises the action will contribute. Again we expect, on the right hand side, 
the exponent to be real negative, i.e., the contribution of a Euclidean path, with i J dt replaced 
by — J dr. This is also reasonable since we know the usual kinetic energy of the particle has to be 
replaced by a negative contribution when the trajectory is under the barrier. 

The summary of this discussion is that actually we should be looking only for the imaginary 
part of the contributions on both the sides of the formula above. If we find the path which 
minimises the Euclidean action, then in terms of that, to leading order, and in the semi-classical 
limit, we have the tunneling formula 

T = Aexp(-S E ) 

We can also now see the reason for T/2 to be the appropriate time. If the path minimises the 
action it should be as symmetric as possible. Thus we expect a time symmetry T — > — T and this 
explains why the escape point <p2 should occur at T/2. We therefore solve the Euler-Lagrange 
equations derived from the action 

»-/*«{Ks) ,+ 5 |v ' l ' +m 

Now the action will be minimum for the path that has the fewest wiggles, i.e., is mostly monotonia 
We expect the path to start at large negative t at the value </>i and stay at that value as much as 
possible, and monotonically reach </>2 near the origin, and then retrace a symmetric path back to 
4>x as r goes to infinity. Such a path which bounces back has been called "the bounce" . 

To solve for the bounce, our first simplification will be to invoke space-time symmetries, viz., 
we assume that the configuration of fields which will minimise the integral in question will obey 
spatial isotropy. This is same as assuming that the spontaneously formed bubble will be spherically 
symmetric. With 4 Euclidean dimensions, assuming 0(4) symmetry, one solves the equation 

f + V - V'{(f>) = 
r 

One boundary condition is <p(r — ► oo) = <j>\ where we have chosen <f>i = 0. At the origin we could 
place the requirement <fi(r = 0) = <f>2 but it is more important to require <j>'(r = 0) = as is usual 
for spherical coordinates when the solution is expected to be smooth through the origin. Indeed 
we may not even know what the "exit point" 02 after the tunelling will be. The bounce when 
solved for will also reveal it. A typical bounce solution is shown in figure 

6.1.1 Next to leading order 

According to above discussion, the tunneling rate is given by a WKB type formula Y = A exp (—Se)- 
With 

SE(4>bounce) determined by extremising the Euclidean action, the exponential is the most 
important factor in this formula. The front factor A arises from integration over the small fluc- 
tuations around the stationary point (^bounce- This is a Gaussian integration and results in a 
determinant. There are several subtleties which arise. The answer is that we need to remove 
zero mode(s) of the fluctuation operator, and normalize with respect to the determinant of the 
fluctuations in the absence of the bounce. Thus with prime on "det" denoting removal of zero 
mode(s), 

( S E (<P) \ 2 ( det' [-D E + V" \ 
~ V 2tt ) \ det [-n E + V" (0)] J 
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Figure 5: Tunneling rate is determined by the minimum of the Euclidean action cj)bounce{r) obeying 
appropriate symmetry and boundary conditions. 

6.1.2 Thermal bounce 

We can now address the tunneling problem in a thermal ensemble, i.e., at T ^ 0. Recall our ob- 
servation during the discussion of thermal effective potential, viz., the analogy between Euclidean 
path integral and the trace weighed by density matrix in thermal partition function. By the same 
arguments we can show that we must look for a bounce solution periodic in imaginary time with 
period 1/(3. There will be an infinite number of such bounces, but the one to dominate will have 
least number of extrcma. Then we obtain the rate formula 



To obtain the euclidian action in this case, we solve the equations of motion obtained by varying 
this action and solve them subject to 0(3) symmetry of spatial directions. This is called the 
"bounce" solution relevant to thermal transitions and the corresponding value of the action is 
inserted into the rate formula. 

6.2 Applications 

The formalism developed in this section is important for determining the evolution of the Universe 
when the field theory signals a first order phase transition. The rate T can vary greatly due to 
the exponential factor. If the rate is too small, the expansion rate of the Universe may be faster 
and in this case parts of the Universe may never tunnel to the true vacuum. If the state of broken 
symmetry is phenomenologically desirable, the rate should be fast enough, at least faster than the 
expansion rate of the Universe at the time of Big Bang Nucleosynthesis (BBN). The end of a first 
order phase transition dumps a certain amount of entropy into the Universe, similar to latent heat 
in usual substances. Since the history of the Universe after the BBN is fairly precisely known any 
unusual phenomenon especially one that may disturb the baryon to entropy ratio should occur 
well before the BBN and not alter the required value of the ratio. 

Such considerations place constraints on the parameters of the scalar field theory undergoing 
the phase transition. The first proposal of inflationary Universe required a fairly specific range for 
the value of the rate, slow enough for sufficient inflation to occur, but still fast enough that the 
present day Universe would be in the true vacuum. 
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In some theories the false vacuum is phenomenologically the desirable one. In many supersym- 
metric models, the desirable state actually turns out to be metastable, whereas the true ground 
state has undesirable properties such as spontaneous breaking of the QCD colour symmetry. 
Several models of supersymmetry breaking arising in a hidden sector and communicated to the 
observed sector by messengers, end up with unphysical ground states. In such cases one invokes 
the possibility that the parameters of the theory make the tunneling rate much smaller than the 
expansion rate of the universe till the present epoch. If a volume of the size M~ 3 determined by 
the energy scale M of the high energy theory is to not undergo a transition within the typical 
expansion time scale of the Universe, then 

r < M 3 H 

where H is the present value of the Hubble parameter, i.e., the experimentally observed Hubble 
constant. 

As another application, in Standard Model, the Higgs boson has a self interaction potential 
of the type discussed above. The exact form of the potential, namely, which local minimum is 
energetically favored and what can be the tunneling rate for going from one vacuum to another is 
determined by the mass parameter and the quartic coupling occurring in the Higgs potential. The 
mass parameter is known from the requirement of spontaneous symmetry breaking to reproduce 
the Weak interaction scale. However, the quartic coupling will be determined only when the 
collider experiments will determine its mass. Since this value is as yet unknown, we can use 
cosmology to put a bound on its possible values. 

It can be shown that if the Higgs boson is very light, then there is a danger for the Universe to 
be trapped in an unphysical vacuum. This puts a lower bound on the Higgs boson mass at about 
10 GeV. This is known as the Linde- Weinberg bound. 

7 Inflationary Universe 

It is remarkable that the Friedmann- Robertson- Walker model is so successful a description of the 
observed Universe. At first this seems a resounding triumph of General Relativity. It is true 
that the dynamics all the way back to Big Bang Nucleosynthesis (BBN) is successfully described. 
However one begins to notice certain peculiarities of the initial conditions. First of all, the Big 
Bang itself presents a problem to classical physics, being a singularity of spacetime. But we expect 
this to be solved by a successful theory of Quantum Gravity. But now we will show that even the 
conditions existing after the Big Bang and well within the realm of classical General Relativity 
pose puzzles and demand a search for new dynamics or newer laws of physics. 

7.1 Fine tuned initial conditions 

Any physical entity such as a planet or a star or a galaxy has an associated legth or mass scale. 
Very often this is set by accidental initial conditions. However not all accidentally possible scales 
would be tolerated by the dynamics holding together the object. Some initial conditions will 
lead to unstable configurations. Question : what scales would be permitted by a universe driven 
according to General Relativity? The coupling constant of Gravity is dimensionful, G^ 1 = M Pl <~ 
(10 19 GeU) 2 . The available physical quantity is energy density. We propose a naive possibility 
that the system size or scale be decided by this energy density re-expressed in the units of Mpi. 
We then find that the Universe has far too small an energy density 10~ 6 (eU) 4 ~ 10~ 66 (Mpi) 4 , 
and far too big a size, Giga parsec compared to 10~ 20 fermi. It is interesting that Gravity permits 
a viable solution with such variance from its intrinsic scales. 

Further, we find that above mentioned ratio evolves with time because the scale factor grows 
and density keeps reducing. We can ask what is the time scale set by gravity for this variation. 
In the evolution equation for the scale factor R(t) we see that the intrinsic scale is set by the 
Gravitational constant. If the p/(Mpi) A was order unity at some epoch, then there is only one 
independent time scale left, that set by Mpi. We then find that the universe would either have 
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recollapsed or expanded precisely within the time set by the Planck scale, 10 sec. The fact 
that the Universe seems to be hovering between collapse and rapid expansion even after 14 billion 
years, requires that we must start with extremely fine tuned initial conditions. The fine tuning 
has to be to the extent of one part in 10 66 because we started with the value of the ratio close 
order unity. 

It is such uncanny fine tuning that is the motivation for proposing an "inflationary" event in 
the early Universe, a phase of unusually rapid expansion. Such an event reconditions the ratios we 
discussed above. Their large apparant values then arise from dynamics rather than initial values. 
All physical realizations of this proposal have consisted of admitting dynamics other than Gravity 
to intervene for the purpose of significantly reconditioning these ratios. 

It is fully likely that the answer to the puzzles to be described in the following is buried in 
the Planck era itself. Any data which can throw light on such a mechanism would also provide 
a valuable window into Planck scale physics. But there are viable candidates within the known 
physical principles of Relativistic Field Theory, a possibility which if true would reduce our intrinsic 
ignorance of the physical world, and in turn lead to prediction of newer forces. 

7.2 Horizon problem 

This problem arises from the fact that our Universe had a finite past rather than an indefinitely long 
past. A finite past gives rise, at any given time, to a definite physical size over which information 
could have travelled upto that time. This physical scale is called the "particle horizon" . At 
present epoch we find unusually precise correlation in the physical conditions across many particle 
horizons, i.e., over regions of space that had no reason to be in causal contact with each other. 
It is to be noted that the world could well have emerged all highly correlated from the Planck 
era. But as explained above, we work in the spirit of exploring newer dynamics within the known 
principles. 

With this preamble, the paradox presented by the current observations is as follows. We know 
that the observed Cosmic Microwave Background Radiation (CMBR) originated at the time of 
decoupling of photons from the partially ionized Hydrogen. The temperature of this decoupling, 
as we estimated in sec. 14.41 is 1200K, while today it is close to IK. From the ratio of temperatures, 
and assuming a matter dominated Universe, 

To 1_ = RjUec) = (td\ m 

T dec ~ 1200 R(t ) \t J 
Therefore, td ~ 2 x 10 5 h^ 1 years. 

Now consider the size of the particle horizon at these two epochs, ie., the size of the region 
over which communication using light signals could have occurred since the Big Bang. 

a (t ) / — — « M « 6000 h' 1 Mpc 

where the contribution of the lower limit is ignorable. Similarly at tdec, horizon ~ Stdec ~ 
0.168 h Mpc, using the time-temperature relation appropriate to the radiation dominated era. 
Then the angle subtended to us today by a causally connected region of the decoupling epoch is 

168 

6 d = ~ 0.03 rad or 2 deg 

6000 5 

This means that we are viewing today Tg^p ~ 14, 000 causally unconnected horizon patches, and 
yet they show remarkable homogeneity. 

7.3 Oldness-Flatness problem 

An independent puzzle arises due to the fact that the curvature of the three dimensional space is 
allowed to be non-zero in General Relativity. Let us rewrite the evolution equation by dividing 
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out by H(t) 2 = (R(t)/R(t)) 2 

1 k 

H{ty R{ty = m ~ 1 

where il(t) = 8irGp(t)/3H(t) 2 which can be thought of as the energy density at time t re-expressed 
in the units of G and H(t). At the present epoch to in the Universe, the observational evidence 
suggests the right hand side (RHS) of above equation is ±0.02. This suggests the simple possibility 
that the value of k is actually zero. Let us first assume that it is non-zero and assume a power 
law expansion R(t) = Rot n , with n < 1 as is true for radiation dominated and matter dominated 
cases. After dividing the previous equation on both sides by the corresponding quantities at 
present epoch to, we find 

— = 50 x (O(i) - 1) 
toy 

where we have used the current value of the right hand side (RHS), 0.02. Now the current value 
t is w 5 x 10 17 second, while at the time of Big Bang Nucleosynthesis (BBN) it was only about 
100 seconds old. Using n = 1/2 for the sake of argument, we find LHS« 10~ 7 , which means that 
correspondingly, on the RHS fl must be tuned to unity to one part in 10 8 . If further, we compare 
to earlier epochs such as the QCD phase transition or the electroweak epoch, we need higher and 
higher fine tuning to achieve the 0.02 accuracy at present epoch. We thus see that the initial 
conditions have to be fine tuned so that we arrive at the Universe we see today. Equivalcntly, 
since the problem is connected to the large ratio of time scales on the LHS, we may wonder why 
the Universe has lived so long. This may be called the "oldness" problem. 

In case the k is zero, then that would be a miraculous fine tuning in itself. The discrete values 
0,±1 arise only after scaling the curvature by a fiducial length-squared. The natural values for 
the curvature pass smoothly from negative to positive and zero is only a special point. The tuning 
of the value to zero earns this puzzle the name "flatness" problem. 

We can also restate the problem as there being too much entropy in the present Universe. 
Taking the entropy density of the CMB radiation at 2.7K and multiplying by the size of the 
horizon as set by approximately Hq 1 <~ (l/3)to ~ 3 x 10 9 yr. we get the entropy to be 10 86 ( 
check this!). The aim of inflationary cosmology is to explain this enormous entropy production as 
a result of an unusual phase transition in the early Universe. 



7.4 Density perturbations 

A last important question to be answered by Cosmology is the origin of the galaxies, in turn of life 
itself and ourselves. If the Universe was perfectly homogeneous and isotropic, no galaxies could 
form. Current observations of the distribution of several million galaxies and quasars suggests 
that the distribution of these inhomogeneities again shows a pattern. To understand the pattern 
one studies the perturbation in the density, Sp(x, D) = p(x, D) — p, where p is the average value 
and p(x, D) is the density determined in the neighborhood of point x by averaging over a region 
of size D. It is found that the density fluctuations do not depend on the scale of averaging D. It 
is a challenge for any proposal that purports to explain the extreme homogeneity and isotropy of 
the Universe to also explain the amplitude and distribution of these perturbations. 



7.5 Inflation 

The inflationary universe idea was proposed by A. Guth to address these issues by relying on the 
dynamics of a phase transition. The horizon problem can be addressed if there existed an epoch 
in the universe when the particle horizon was growing faster than the Hubble horizon. Later this 
phase ends and we return to radiation and matter dominated Universe. 

If the Universe was purely radiation dominated during its entire early history particle horizon 
could not have grown faster than Hubble horizon. However, if there is an unusual equation of 
state obeyed by the source terms of Einstein's equations then this is possible. Our study of first 
order phase transitions suggests a possible scenario. We have seen that there is a possibility for 
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the system ( the universe) to be trapped in a false vacuum. Exit from such a vacuum occurs by 
quantum tunneling. If this tunneling rate is very small, the vacuum energy of the false vacuum 
will dominate the energy density of the Universe. Vacuum energy of a scalar field has just the 
right property to ensure rapid growth of particle horizon, keeping the hubble horizon a constant. 
If inflation occurs, the flatness problem also gets automatically addressed. 
Consider the energy momentum tensor of a real scalar field 

Tj? = d^d^-S^C 



= d^d^ - Si QdA0d A </> - V^) 

If the field is trapped in a false vacuum, its expectation value is homogeneous over all space and is 
also constant in time. This means all the derivative terms vanish and Tjf —* VqS^ where Vq is the 
value of the potential V in the false vacuum, called the vacuum energy density. Now for an isotropic 
and homogeneous fluid, the energy- momentum tensor assumes a special form, diag(p, — p, — p, —p). 
That is, the off diagonal terms are zero, the three space dimensions are equivalent, and the non- 
zero entries have the interpretation of being the usual quantities p the energy density and p the 
pressure. Thus comparing this form with that assumed by the scalar field in a false vacuum, we 
sec that the expectation value of the scalar field behaves like a fluid obeying the unusual equation 
of state p — —p. Now the Friedmann equation becomes 




whose solution is R(t) = R {U) exp (H (t - U)), with H 2 = (87r/3)F and U is some initial time. 



7.6 Resolution of Problems 
7.6.1 Horizon and flatness problems 

During the inflationary epoch we assume an exponential expansion and a constant value of Hubble 
parameter H (defined without special subscript or superscript since Hq is reserved for the current 
value of H). Now consider the particle horizon, or equivalently, the luminosity distance at any 
epoch t 

*—"/,.£ "iff"™ 

Thus the distance over which causal effects could be exchanged is exponentially larger than the 
simple estimate Mdec we used while discussing the horizon problem. Suppose inflation lasted for 
a duration (t — U) = r. We can estimate what value Ht we need during inflation so that we are 
not seeing a large number of primordial horizon volumes but only about one. Suppose inflation 
ended leaving the Universe at a temperature T r ( subscript "r" signifies reheat as explained in 
next subsection). Then the dn of above equation rescaled to today assuming radiation dominated 
Universal should give current inverse horizon Hq 1 . Thus 

1 u T \ fT r \ 1 



— e x , 

H ) \T J H 

In above formula let us estimate H/Hq again assuming radiation dominated evolution since 
the T r till now, 



FRW 



TrV / 10 14 GeV V 
YJ ~ I 10- 4 eV ) 



do 27 ) 



4 The Universe has of course not been radiation dominated through out. But replacing later history of the 
Universe by matter dominated evolution will not significantly alter these estimates since the expansion is changed 
by a small change in the power law, t 2 / 3 instead of t 1 / 2 . Assumption of radiation dominated expansion allows 
relating temperatures at two different epochs, as a good approximation. 
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Figure 6: Comparison of simple FRW cosmology with inflationary cosmology. Figure courtesy 
Kolb and Turner. 
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recall that H 2 oc T A and we inserted a Grand Unification scale 10 14 GeV as a possibility for T r . 
Then e Hr has to be 10 27 . 

We need an improvement upon this estimate. As we shall later see, there are theoretical 
reasons to believe that T r < 10 9 GeV. In this case it is also assumed that there is additional ( 
non-inflationary or mildly inflationary) stretching by a factor 10 9 . In this case we delink the 
reheat temperature T r from the vacuum energy density causing inflation, and assume the latter 
to continue to be at GUT scale. In this case, 

Ht. ( H \ „ (to \ N , i n-9 . . 1 0I8 



where we have used estimate of the vacuum energy density (0.03eV) 4 as derived from the directly 
observed value of current Hubble constant. In the literature one often sees this estimate made 
with GUT scale taken to be 10 16 GeV so that the required value of e Hr > 10 22 ~ e 55 . 

Let us see the requirement to solve the flatness problem. We need to show that the term \k\/R 2 
becomes insignificant regardless of its value before inflation. For naturalness we assume it to be 
comparable to H 2 as expected from Friedmann equation. Now including scaling from inflation, 
the stretching by factor 10 9 as introduced above, and subsequently during the era after reheating, 
we find 



e>2 ; i d2 i 

v Jt /o \ 1L / prc-mt 

The left hand side is (1 — Q)H 2 which is insignificant. The right hand side is the same small factor 
we estimated above, squared. So this means the left hand side is reduced to a value 10 -36 or 10 -44 
depending on the value of Grand Unification scale we take. 



7.6.2 Avoidance of unwanted relics 

A byproduct of Inflation is that it would also explain absence of exotic relics from the early 
Universe. Typically a grand unified theory permits occurance of topological defects such as cosmic 
strings or monopoles. They signifiy unusual local vacua of spontaneously broken gauge theories 
which cannot evolve by unitary quantum mechanical processes to the simple vacuum. Unlike 
normal heavy particle states, therefore, such defects cannot decay. 

The natural abundance of such relics can be calculated by understanding the dynamics of their 
formation. Typically these events are the phase transitions characterised by specific tempratures. 
If the naturally suggested adundances of these objects really occured, they would quickly dominate 
the energy density of the Universe, with possible exception of cosmic strings. This would be 
completely contradictary to the observations. On the other hand if the scale of Inflation was 
below the temperature of such phase transitions, the density of topological objects formed would 
be diluted by the large factor by which volumes expand during inflation. 

Another class of exotic relics are the so called moduli fields, scalar excitations arising in su- 
persymmetric theories and String Theory. They are generic because of the powerful symmtry 
restrictions on potential energy functions in such theories. Inflation provides a solution for some 
class of models for this case also. 



7.6.3 Resolution for density perturbations 

Finally the density perturbations are neatly explained by inflation. The scalar field is assumed to 
be in a semi-classical state. However quantum fluctuations do exist and these should in principle 
be observable. Inflationary era is characterized by Hubble parameter H remaining a constant 
while the scale factor grows exponentially. Thus the wavelengths of various Fourier components 
of the fluctuations are growing rapidly, leaving the horizon. 

We shall take up in greater detail the theory of small perturbations in an expanding Universe 
in section [SJ There we show that the amplitudes of the Fourier modes of these perturbations 
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remain frozen at the value with which they left the Hubble horizon. Eventually when the Universe 
becomes radiation dominated and subsequently matter dominated, Hubble horizon H^ 1 begins to 
grow faster than the scale factor and the wavelengths of the modes begin to become smaller than 
Hubble horizon. This is the same reason as the solution of the horizon problem wherein apparently 
uncorrelated regions of distant space now seem to be correlated. They all emerged from the same 
causally connected region and got pushed out of the horizon during the inflation epoch. 

The result of this evolution of the fluctuations is that when they re-enter the horizon they all 
have the same amplitude. These fluctuations then influence the rest of the radiation and matter 
causing fluctuations of similar magnitude in them. This is the explanation for the scale invariant 
matter density perturbations represented by distribution of galaxies. We shall take up the details 
of fluctuations in the next section. 

7.7 Inflaton Dynamics 

Inflation is a paradigm, a broad framework of expectations rather than a specific theory. The 
expectations can be shown to be fulfilled if a scalar field dubbed "inflaton" obeying appropriate 
properties exists. If inflation is implemented by such a scalar field, we need to make definite 
requirements on its evolution. We assume that its evolution leads the Universe through the 
following three phases 

• Inflationary phase 

• Coherent oscillation phase 

• Decay and re-heating phase ("re-heat" only for low field scenario) 

Of these three phases, the inflationary phase addresses the broadest requirements discussed 
in previous subsection. To obtain exponential expansion we need constant vacuum energy. This 
means that the field phi has a value where V(^) ^ and also that <p continues to remain at such a 
value. The simplest such possibility is a false vacuum, a local minimum of the effective potential 
which is not a global minimum. But this is problematic because this kind of state can be far too 
stable and the exit from it keeping the Universe permanently inflating. Two possibilities which 
are strong candidates due to phenomenological reasons are the so called High Field (older name 
Chaotic Inflation) or the Low Field ( older name New Inflation) scenarios of inflation. These are 
shown in figure [7J In the High Field case the initial value of the field is close to Planck scale and 
no clear barrier separating it from the low energy true minimum. However its dynamics governed 
by Planck scale effects is "chaotic" and keeps it at a very high energy for a long time. In the 
Low Field case, the dynamics is usual field theory but the effective potential function has a long 
plateau of very small slope. Assuming the initial value of the field at the top of the plateau, this 
allows the field to sustain a position of large vacuum energy for a long time. In both cases, the 
field eventually moves towards the low energy true minimum, via the next two phases. A third 
possibility which is appealing for supersymmetric unified theories is called the Hybrid scenario. 
It involves two fields, one which keeps the Universe initially at a high field value, and the second 
field which becomes more dominant at a later stage, causing a rapid roll down and exit from the 
high energy plateau. It will not be possible for us to discuss these scenarios in any detail in these 
notes. 

Next, it is easier to explain the third listed phase, since it is necessary that the end point of 
inflationary expansion is a hot Universe. Big Bang Nucleosynthesis is very successful in explaining 
the natural abundance of elements. We need to assume a hot Universe of at least a few MeV 
temperature for BBN to remain viable. 

Finally, the intermediate phase dominated by coherent^ oscillations of the scalar field is almost 
certain to occur as inflation ends. This is because during inflation the field is already in a coherent 
state, one in which it has homogeneous ( position independent ) value. This phase ends with 
creation of quanta. In a large class of models this phase may be of no particular interest. But 

5 i.e., a special state in which it is possible to treat the field in the leading order as if it were a classical field. 
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Figure 7: Sketches of stipulated effective potentials with the corresponding initial value of the 
inflaton field indicated by a small circle ( a "ball" ready to "roll down" the shown profiles). Left 
panel shows the High Field scenario, the right panel shows the Low Field scenario 

in special cases when the coupling of the inflaton to other matter is tuned to certain values, it 
can lead to a variety of interesting effects which can have observable consequences. One such 
possibility is a long duration of coherent oscillations, which can be shown to mimic a Universe 
filled with pressureless dust. Alternatively there can be particles with special values of mass which 
can be shown to be produced copiously during the coherent oscillation epoch. Such effects are 
called "preheating" because heating of the Universe is achieved not directly through the inflaton 
but by other particles which have efficiently carried away the energy of the inflaton. Such phases 
can enhance the expansion achieved during the inflationary phase, the additional f 9 factor used 
in estimates in previous subsection. 

7.7.1 Predictions and observables 

For a quantitative study we consider a scalar field (f> with lagrangian 

S[4>] = J d^x^—g - V{<j>) + Vnmtter) 

To discuss the inflationary phase we do not need £0_ m atter containing the detail of the coupling 
of the scalar field to the rest of matter. Varying this action after putting the Fricdmann metric 
for gives an equation of motion for <j>, 

4> + 3H<p + V'((f)) = 

The inflationary phase is characterized by a homogeneous value of <j> and a very slow time evolution 
so that there is domination by vacuum energy. Mathematically one demands 

0< 3H<j> 

in other words we assume that the time scale of variation of the field <f> encoded in 0/ " ■ <fi is small 
compared to the time scale H^ 1 of the expansion of the Universe. Further, the assumption that 
the time scale of variation of 4> is very small amounts to assuming that in the action we must have 

2 « U 

Thus while obtaining the Euler-Lagrange equations the V term continues to be important, but 
the higher time derivative of <j> can be dropped. So the evolution equation reduces to 

3H4> = -V'{(f>) 
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To ensure consistency of the above two conditions, divide the simplified evolution equation by 
3H and take a time derivative. Recall that H 2 — (8tt/3)GV (</)). Thus H can be implicitly 
differentiated with respect to time as shown in the equation below 

We now write the equation in three equivalent forms 

4> _ V{<t>) , V 



3Hj> 9H 2 9m 



H' 



V"(d>) + V'(f)GV> 



9xfGVb 9x(fG)V 2 



V" 1 1 (V 

+ - 



3 x 8ttGV 2 3 x 8nG \ V 

where use has been made of the preceding assumptions as also the relation 2HH' — ^-GV. 
Inspecting the above equations we define three parameters 



_ Mp 2 ( V'(<f>) y Mp 2 V" (4>) 

and 



£ ~ 16tt ^ V{<P) ) ' V ~ 8tt V{4>) 



H<f> 

The requirements of the inflationary phase viz., large vacuum energy and a vary slow roll 
towards the true minimum mean that 

e < 1, | f] |< 1 and ^ - 0(e,r/) 

We use these in the criterion that if either of these quantities becomes large, inflationary phase 
ends. These have come to be called the "slow roll" parameters characterizing inflation. 

For several decades inflation remained a theoretical paradigm. However with precision cos- 
mological experiments such as the Hubble Space Telescope (HST) and the Wilkinson Microwave 
Anisotropy Probe (WMAP) yielding valuable data we face the exciting prospects of verifying and 
refining the paradigm, and also deducing the details of the dynamics of the inflaton field. It is pos- 
sible to set up a relationship between slow roll parameters introduced above and the temperature 
fluctuation data of the microwave background radiation. Similarly the large scale galaxy surveys 
such as 2dF GRS ( 2 degree Field Galaxy Redshift Survey) and 6dF GRS provide detailed data on 
distribution of galaxies which can be counter checked over a certain range of wavenumbers against 
the fluctuations as observed in WMAP and also against the specific dynamics of the inflaton. 

Finally let us work out a simple example of how we can relate intrinsic properties of the effective 
potential to observable features of inflation. We can for instance compute the number of e-foldings 
N in the course of <\> evolving from initial value (pi to a value 4>f, given the form of the potential. 
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3 H 2 d<j) 
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where we assumed t/ <~ H 1 (</>/). Now shift the global minimum of the effective potential to be 

at 4> = so 4>i » <pf 




V{<t>) 



d<f> 



V'{4>) 



This formula can be used to relate the dominant power law in the effective potential with the 
number of e- foldings. For V{<p) — A</> M , 



8 Density perturbations and galaxy formation 

An outstanding problem facing FRW cosmological models is formation of galaxies. If the Uni- 
verse emerged from the Planck era perfectly homogeneous and isotropic, how did the primordial 
clumping of neutral Hydrogen occur? Without such clumping formation of galaxies and in turn 
stars would be impossible. 

The related observational facts are also challenging. The fluctuations in the average density 
have resulted in a distribution of galaxies and clusters of galaxies. What is remarkable is that 
these fluctuations exist at all observable scales. Further, the observed fluctuations seem to have 
originated from seed fluctuations which were of the same magnitude, approximately one part in 
10 5 , independent of the scale at which we study the fluctuations. This statement of scale invariance 
began as a hypothesis, known as the Harrison-Zel'dovich spectrum but has been remarkably close 
to the extensive experimental evidence accumulated over the last fifty years. The main sources 
of current data are Sloan Digital Sky Survey (SDSS), Two degree field Galaxy Redshift Survey 
(2dF GRS), quasar redshift surveys, "Lyman alpha forest" data etc, collectively called Large Scale 
Structure (LSS) data. 

Here we shall present a brief overview of the formalism used for studying fluctuations. We 
shall also show that the inflationary Universe is in principle a solution, providing scale invariant 
fluctuations. The magnitude of resulting fluctuations however is too large unless we fine tune a 
parameter to required value. 

8.1 Jeans Analysis for Adiabatic Perturbations 

First we study the evolution of perturbations in a non-relativistic fluid. We study the continuity 
equation, the force equation and the equation for gravitational potential assuming that the state 
of the fluid provides a solution a solution to these. We then work out the equations satisfied by 
the perturbations. Fourier analysing the perturbations, it is found that modes with wavelengths 
larger than a critical value of the wavelength the perturbations are not stable. 

The relevant equations for the mass density p(x, t), velocity field v(x, t) and the newtonian 
gravitational potential 0(x, t) are 
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We now split the quantities into average values p, p and space-time dependent perturbations p\ , 
Pi 

p(x, t)=p + p 1 (x, t) p(x, t)=p + Pl (x, t) 

and similarly for velocity. However, for a homogeneous fluid, average velocity is zero so v and vi 
are the same. It is reasonable to assume that there is no spatial variation in equation of state. 
This means the speed of sound is given by 



v 2 



dp\ pi 
•Wadiabatic Pi 



Thus the equations satisfied by the fluctuations are, 

+ ^Vpi + V<Ai = 
at p 

V 2 0i = 4ttGpi 
From these coupled equations we obtain a wave equation for pi, 

^- v ^ Pl =4nGp Pl 

whose solution is 

pi(r,t) = Ae (iwt " ik - r) 
withw 2 = v 2 s k 2 - 4irGp 
The expression for u> suggests the definition of a critical wavenumber, Jeans wavenumber, kj = 

1 /2 

(AnGp/vj) ' . For k < fe / we get exponential growth, i.e., instability. We can understand 
this result by associating with a wavelength A a hydrodynamic timescale T^yd ~ X/v s ~ l/kv s . 
This is the timescale during which pressure differences will be communicated by perturbations 
with wavelength A. Next we associate the timescale Tg ra v = (^Gp ) ' with the gravitational 
influences. The result above says that for the wavelengths for which hydrodynamic response is 
slower to propagate than gravitational influences, the latter win and cause a gravitational collapse. 
We can now define Jeans mass 

„ 47r/7r\ 3 tt 5 / 2 v s 3 
Mj = — — Po 



3 V^/ 6 G 3/2 p l/2 

We can deduce that a homogeneous mass bigger than this value is susceptible to gravitational 
collapse. 

8.2 Jeans Analysis in Expanding Universe 

We need to extend the above treatment to the case of expanding Universe. Firstly we must let 
the mean density and pressure be time dependent due to change in cosmological scale factor. For 
instance, the mean density in a matter dominated universe will scale as p(t) = p(U) (R(U / R(t)) 3 
where U denotes some "initial" reference time. 

From now on we shall not follow the evolution of the other quantities but focus on the energy 
density. We introduce the dimcnsionlcss quantity 5 = pi/p from which the obvious R(t) depen- 
dence gets scaled out. Further, using the FRW metric in the comoving form, we treat r to be 
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dimcnsionless and t and R(t) to have dimensions of length ( equivalently, time). Introduce the 
Fourier transform 

5(x ' f)= /p ik( ' )eHk ' r) 
It can be shown that these Fourier modes obey the equation 

OP, /v 2 k 2 \ 

5 k + -fS k + ( M- - 4nGp(t)) 4 = 

We see that the modified Jeans wave number defined by 

k 2 ee AnGp(t)R(t) 2 /^ 

plays a crucial role, in that at any given epoch t, wavelengths shorter than ~ are oscillatory 
and hence stable. 

We shall now study the fate of the long wavelengths, the ones significant on cosmological scales. 
We shall see that the unstable modes grow in time, but instead of exponential growth they can 
have power law growth. We have assumed k — > 0. Further, using FRW equation we can replace 
the Gp{t) term by {3/2){R/R) 2 in a spatially flat universe. Then for matter dominated universe 
with R oc t 2 / 3 , we get 

4 • 2 

S + ttS - ^5 = 

The solutions are 

6+(t) - S (U)(- 



6-(t) = 5-{U)(£) 

Studying the examples of other power law expansions of scale factor R, one may conclude 
that the expansion of the universe keeps pulling apart the infalling matter and slows down the 
growth of Jeans instability. However In the case of de Sitter universe it is found that exponential 
instability persists. In this case the long wavelength modes obey the equation 

5 + 2H5- '^H 2 S = 

so that substituting S <~ e at we find 

a 2 + 2Ha - hi 2 = 

This has the roots a± = —H ± y/ H 2 + 3/2. One root is negative definite signifying decaying 
exponential though one positive root persists, a + = \J H 2 + 3/2 — H. 



8.2.1 Fate of the super-horizon modes 

Let us now return to the idea of inflationary universe as the source of primordial perturbations. The 
basic hypothesis is that the quantum mechanics of the "inflation" scalar field causes fluctuations 
in its expectation value. These then manifest as perturbations in the classical quantity, the energy 
density. This assumption is expressed as 



Pi = Sp = 




where 

5<t>EE({4>-m 2 ) 
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and the expectation values are computed in an appropriately chosen state. This choice is not always 
easy. In general this is a static, translation invariant state with similar properties shared by the 
vaccum expectation values. In an expanding Universe the corresponding symmetries available are 
those of the space-time metric, namely the FRW metric. Since the time translation symmetry 
is lost, there are several conceptual issues. Fortunately the de Sitter metric has a sufficiently 
large group of symmetries permitting a fairly unique choice of the vaccum. Inflationary universe 
resembles the de Sitter solution over a substantial length of time so that we can adopt the answers 
obtained for the de Sitter case. 

Decompose 5(j> into Fourier modes with the same conventions as in the preceding section. 
During the inflationary phase, the expectation value of <j> remains approximately constant. Hence 
by appropriate shifting of the field, dynamics of 5(j) and <j> are the same. The equations of motion 
for the modes of S(j) are then 



Then for "super-horizon" fluctuations with wavenumbers satisfying k RH, we can ignore the 
last term and the non-decaying solution is S(f>k = constant. This is a crude argument to justify 
that the fluctuations at this scale become constant in amplitude. But constant amplitude would 
mean vanishing time derivatives, so the third term can't be smaller than the first two. In order to 
consistently ignore the third term relative to the second, we need to additionally assume that if the 
time scale of variation of 54>k is r, then 1/r > k/R in addition to k <C RH. The two inequalities 
together imply 1/r S> H. Thus the fluctuations are constant in the sense that the time scale of 
their variation is much less than the natural time scale of the geometric background, if -1 . 

The above statement can be made more precise using Quantum Field Theory in curved space- 
time, where it can be shown that in de Sitter universe, for a massless scalar field, the fluctuations 
in <fi, after appropriate cut-off procedure, are given by 



The assumption in this calculation is that the state chosen is de Sitter invariant. Thus variations 
that do not respect this invariance don't constribute and we get the expected constant result. 
Further, since no new scales are introduced through the choice of the state, the only dimension- 
ful quantity available in the problem is H, which effectively determines the magnitude of the 
fluctuations. 

We can now see the qualitative features which make inflation so appealing for generating scale 
invariant perturbations. For all normal kind of states of matter and energy, the scale factor grows 
as power law t s , with s < 1. Thus H^ 1 grows as t and all wavelengths X p h y scale as (X p h y /R (ti)) x 
R(t) and hence keep falling inside the horizon. On the other hand, in the inflationary phase, H^ 1 
= constant and wavelengths grow exponentially, X p h y {t) oc Aje^* - **) with some reference initial 
time tj. Later, after inflation ends, H^ 1 begins to grow faster than other length scales and steadily 
catches up with fluctuations of increasing values of wavelengths. Between leaving H~ x and re- 
entering H{£)~ x later, the amplitude of the fluctuations remains frozen by arguments of preceding 
para. Note that for other power law expansions t s with s > 1 also, the wavelengths grow faster 
than the Hubble horizon. However, the amplitude of these wavelengths will not remain constant 
and will not reproduce the scale invariant spectrum after re-entering the horizon. Some of these 
arguments will become clearer in the following subsection. 

8.2.2 Connection to density and temperature perturbations 

The problem of galaxy formation is to predict the observed pattern of clumping of luminous 
matter. Also, since matter and radiation were in equilibrium upto decoupling, the fluctuations in 
matter have also to be reflected in the fluctuations in the temperature of the CMBR. 

The fluctuation we are referring to are in the spatial distribution. These are mathematically 
characterized by the auto-correlation function 



S4> k +3H5<j) k + k 2 ^ =0 




(Ap) 2 (r)E(Mx)Mx + r)) 
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where on the right hand side, an averaging process is understood. For a homogeneous medium, 
the locations x are all equivalent and this dependence drops out at the end of averaging process. 
Introducing the Fourier transform <5k, we can show that 



Now the rms value Ap rms at a given point is the square-root of this auto-correlation function for 
r = 0. Accordingly, taking the limit kr — > in the above expression, we get 



V represents the variation in (Ap rms ) 2 with variation in lnfc. The aim of experiments is to 
determine the quantity V. 

Linear perturbation theory used above is valid only for small fluctuations. Once a fluctuation 
grows in magnitude it begins to be controlled by non-linear effects. We can estimate the intrinsic 
scale upto which the mass of a typical galaxy could have been in the linear regime. Using the 
value of the p to be the present abundance of non-relativistic matter ( ss 10 _29 g/cc) and bringing 
out a factor of 10 11 solar masses, the mass in a sphere of diameter of wavelength A is 



Assuming 10 12 solar masses per galaxy, this gives the size A to be 1.9 Mpc, far greater than the 
actual galactic size 30 kpc. The A found here represents the size this mass perturbation would 
have had today had it not entered the non-linear regime. 

Present data are not adequate to determine the spectrum V over all scales. However too large a 
magnitude of fluctuations at horizon scale would have been imprited on temparature fluctuations of 
CMB, which it is not. Likewise large fluctuations at smaller scales could have seeded gravitational 
collapse and given rise to a large number of primordial black holes, which also does not seem 
to be the case. Hence the spectrum must not be varying too greatly over the entire range of 
wavenumbers. It is customary to assume the spectrum of |<5k| 2 to not involve any special scale, 
which means it must be a power law k n . Further, a fluctuation of physical scale A contains mass 
M <~ A 3 <~ fc~ 3 . Hence the spectrum V ~ M- 1 -"/ 3 . Now if we make the hypothesis that the 
spectrum of perturbations seems to be independent of the scale at which we observe it, we expect 
V <~ a constant, i.e. independent of M. For this to be true, n must be —3. 

In the analyses of WMAP data of CMB it is customary to normalize the rms perturbation 
spectrum by its observed value at (8/h) <~ UMpc and denote it as- WMAP reports the best 
fit value to be 0.8. Further, (Ap rms ) 2 is parameterized as fc( _1 + ns ) where the subscript in n s 
signifies scalar perturbations. A very long epoch of perfectly de Sitter inflation would produce 
n s = 1 and a perfectly scale invariant spectrum. Given a specific model of inflation the small 
departures of n s from unity can be calculated as a function of k. This mild dependence of n s on k 
is referred to as " running of the index" of the power law. Current WMAP data ( ie horizon scale 
perturbations) seem to suggest n s = 0.97 but the direct observations of LSS data ( galaxies and 
clusters of galaxies) suggest n s > 1. 

Likewise a formalism exists for relating the temperature fluctuations with the density pertur- 
bations and in turn with the scalar field fluctuations. We shall not go into it here and the reader 
is referred to the references. 

8.3 Density Fluctuations from Inflation 

We now show how the inflation paradigm along with the knowledge of the form of the scalar 
potential helps us determine the magnitude of the scalar field fluctuations. 





with 



2ir z 
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Suppose we wish to know the fluctuation in a scale of size of our present horizon. According to 
the derivation in previous subsection we need to know the value of the perturbation when it left 
the horizon in the inflationary era. And we need to know the number of e-foldings the inflationary 
universe went through before becoming radiation dominated. It is the latter fact which then 
determines the later epoch when the same scale re-enters the horizon. 

Let us trace a physical scale £q today by keeping track of corresponding co-moving value I. 
We have to consider the evolution in two parts. From the present we can only extrapolate back 
to the time when the current hot phase of the Universe began, i.e. the "re"heatec0 phase. Prior 
to that was the phase of inflaton oscillation and decay. Reheating is assumed to be complete at a 
time td, the decay lifetime of the inflaton or its product particles. 

Thus the size of a scale I can be extrapolated to the epoch tf when inflation ended, (i.e., the 
slowness conditions on the evolution of the scalar field ceased to be valid) by 

'To} (R{t f ) 



tf *\T r )\R(t d ) 

The last ratio can be estimated if we are given the effective potential V((f>) = \(j) v and a formalism 
for the dissipation of the inflaton vacuum energy. We shall not pursue these details here but claim 
that this can be calculated to be 

R(t f ) f tf \^^ (1) 



R(td) \td 



1/2 



for v = 4 (2) 



Now tphys (tf) = H^ 1 e e where Nt is the number of e-foldings between the time the specific 
scale attained the value -ff -1 , (i.e. became comparable to the horizon) and the end of inflation. 
Substituting the current value of the horizon in the above expressions finally gives N H -i « 50 — 60. 

We now trace the magnitude of the perturbation through this exit from horizon followed by 
the re-entry at present epoch. It turns out that Sp/ p is a physically ambiguous quantity to follow 
through such an evolution. This is because choice of a particular time coordinate amounts a choice 
of a gauge in General Relativity. The gauge invariant quantity to focus on has been shown to be 
( = 5p/(p + p). 

We seek the value of the numerator at a late epoch when inflation has ended. The denominator 
at this epoch is determined by p = and the energy density which is dominated by the kinetic term. 
The value of £ at the inflationary epoch is known from preceding arguments about perturbations 
on scales comparable to horizon. Here 

^=^ = ^)| 

where the Sip is estimated from QFT calculation of the rms value. Further, we replace </> 2 by using 
the slow roll condition of inflation, 3Hcf> = —V'(<f>). Thus we equate 



-if- 1 



ei^if- 1 

V'(<f>)H(<t>) 



2lT(f> 2 

9ff 3 (0) 
2n(V'(cf>)) 



-ff- 



-ff- 



6 We remind the reader that it is possible in some inflationary scenarios for the Universe to never have been in 
thermal equilibrium before this stage. Hence the prefix "re" is purely conjectural though conventional. 



41 



Thus in matter dominated era when p = 0, we have recovered 



5 J M*V'(<t>t) 

where we reexpress G = 1/Mp, the squared inverse of the Planck mass in natural units. The 2/5 
factor is acquired during transition from radiation dominated to matter dominated era. 

We thus need the values of V and V at the value fa . We do not really know the latter directly. 
But we can determine it if we know the number of e-foldings between its crossing the horizon and 
the end of inflation. Inverting the relation 



H- 1 (fa) = £ phy (t f )e N ' 



for fa and also using 



h) = j Hdt = j 



N e (fa -» 4> f ) = Hdt= ^d<j) -> ttG(/) 



where the last arrow gives the answer corresponding to the form X<f> for the effective potential and 
is obtained by consistently using the slowness condition and the FRW equation. This gives us the 
number of c-foldings between horizon crossing by this scale and the end of inflation. Therefore, 
trading Ne for fa we get 

3/2 



p ) n 5 " v m p ) 5 * \ n ) 

We have arrived at a remarkable mathematical relationship, expressing the magnitude of per- 
turbations visible in the sky today with the parameters of the effective potential that drove the 
primordial inflation. We can take the fluctuations Sp/p to be as visible in the CMB temperature 
fluctuations, ST/T ~ 6 x 10 . Assuming N — 55 as needed for solving the horizon and flatness 
problems, we find that we need the value of A ~ 6 x 10 14 . This is a tremendous theoretical 
achievement. Unfortunately the required numerical value is unnaturally small and it appears that 
we have to trade the fine tuning required to explain the state of the Universe with a fine tuning 
of a microscopic effective potential. 



9 Relics of the Big Bang 

As the Universe cools reaction rates of various physical processes become slow. When they become 
slower than the expansion rate of the Universe, the entities governed by those reactions no longer 
interact and remain as residual relics. The mathematical description for these events is provided 
by Boltzmann equations which can be used to infer the relative abundance of these relic particles 
which can be in principle observed today. The term relic applies to a wide variety of objects, 
including extended objects like cosmic strings or domain walls but we shall be dealing only with 
particle like relics in these notes. 

9.1 Boltzmann Equations 

Boltzmann equations describe the approach to equilibrium of a system that is close to equilibrium. 
In the context of the early Universe we have two reasons for departure from equilibrium. One is 
that if the reheat temperature after inflation has been T re h then all processes requiring energies 
larger than T are suppressed and such processes play no role in establishing dynamical equilibrium. 
Thus particles that interact via only such processes remain out of equilibrium. We do not have 
much more control on such entities and in any case they most likely got inflated away and will 
not be recreated due to insufficient energy for them to be created. 

The more important source of departure from equilibrium is the fact that the Universe is ex- 
panding. Like in an expanding gas, the temperature systematically falls. The primary assumption 
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here is adiabaticity - i.e. extreme slowness of the rate of change of temperature compared to the 
time scales of the equilibrating processes. However, interesting epochs in the Universe correspond 
to times when the rates of a few specific processes are becoming as small as the expansion rate 
of the Universe. After the epoch is passed same reactions go out of equilibrium and the last 
conditions remain impritned as initial conditions for the rest of the evolution. 

Schematically one can think of the Boltzmann equation as Liouville operator L acting on 
distribution function /, with a driving force provided by a collision operator C. In the absence of 
the collision term we have equilibrium statistical mechanics. L[f] — C[f] The Liouville operator 
which basically describes convection through the phase space can be written as 

t d ^ - F ~ 
L = — +v-W x + — ■ V„ 
at m 

assuming the conjugate momentum has the simple form of velocity times mass. In General Rela- 
tivity this has to be generalized to 



P dx * h*P P dpa 



which simplifies in the FRW case to 



E dt~R^ dE 



The total number density is obtained by integrating the distribution function over all momenta, 

.9 

(2tt) 3 



Thus we obtain 



3 df R g f 3 \p\ 2 8f g f nlf] d 3 p 



(2tt) 3 J 1 dt R (2tt) 3 J 1 E dE (2tt) 3 J w J E 

Exchanging the order of integration and differentiation in the first term and working out the 
second term by doing an integration by parts, we can show that this equation becomes 

—n + 3—n = 1 ' ' 



dt R (2?r) 
Consider a process involving several particle species ip, a, b... 

ip + a + b <-> i + j + ... 

Our interest is usually a specific species which is undergoing an important change. We think of 
the collision operator with species tp as the object of main interest to be 

9 „ r , 1 d 3 W' f d?pi P d 3 p a 4 4 



(2tt) 3 lJi dE^ J (2tt) 3 2^ (27r) 3 2S a 



\Mt.Ufa{l±fa) 

- \M\l_ ft (1 ± f a ) (1 ± U) 

where M. represents a matrix element for the concerned process. There are Bose-Einstein and 
Fcrmi-Dirac distribution functions for the species in the in state. As for the out state, the ± 
signs have to be chosen by knowing the species. Bosons prefer going into an occupied state (recall 
harmonic oscillation relation a^\n) = \/n + l\n) so that an n-tuply occupied state has weightage 
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proportional to n to be occupied). Hence the factors (1 + /), while fermions are forbidden from 
transiting to a state already occupied hence the Pauli suppression factors (1 — /). There is an 
integration over the phase space for each species. 

Let us specialize the formalism further to cosmologically relevant case with isentropic expan- 
sion. Since the entropy density s scales as 1/i? 3 , we can remove the presence of R 3 factors in 
number density of non-relativistic particles by studying the evolution of their ratio with s. Thus 
we define the relative abundance for a given species Y = Then the we can see that 

h$ + 3Hn^, = sY 

Next the time evolution can be trade for temperature evolution in a radiation dominated universe, 
by introducing a variable x. 

x = ^so that t = 0.3.g* ^"T^p" = 0.3g* 1 ^ 2 —^-x 2 = H~ 1 (m)x 2 

Thus we get the equation 

dn^dwa, ...G?7r i ...(27r) 4 |M| 2 (5 4 (p in - p out ) 



dY 

dx 



H(m 



I 



[fa fb---fip fi fj---] 

Consider a species ip which is pair annihilating and going into a lighter species X, tp^ — > XX. 
The assumption is that the X are strongly interacting either directly with each other or with rest 
of the contents, so that they equilibrate quickly and remain in equilibrium. Thus the species to be 
studies carefully is ip. Due to the property of the chemical potential and detailed balance which 
would exist if the tp are also inequilibrium we can relate the equilibrium values 



n x n x =n-friQ- = J 



eq = ^ EQ y 



Note the superscript eq is not necessary in the nx due to it always being in equilibrium. We can 
thus obtain the equation 



The solution of above equations can be simplified by identifying convenient regimes of values of x 
in which approximate analytic forms of Yeq exist 



Yeo(x) = ——r I — I — x 3 ' 2 e~ x x ^> 3 non-relativistic case 
27T 4 V8/ 3* s 



Y EQ (x) = -^-£(3)— = 0.278^£ x < 3 re lativistic case 
27r 4 g* s g* s 

where the effective degeneracy factors g e fj are defined relative to their usual values g by g e ff — 

gboson and geff — 4 5 'f 'ermi- 



9.1.1 Freeze out and subsequent evolution 

We can get further insight into the special case considered above, namely that of a species anni- 
hilating with its anti-particle and also going out of equilibrium. Define 

= n E Q (<taW\) 
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which represents the rate of the reactions, given as a product of the microscopic cross-section a, 
and number density times relative velocity as a measure of the flux. Using this, we can rewrite 
the evolution equation above in the form 

x dY _ Fa 
Y EQ dx H 

This shows that the rate of approach to equilibrium depends on two factors. The second factor is 
the extent of departure from equilibrium, as we may expect even in a laboratory process. The front 
factor r a/ H represents the competition between the annihilation rate ( temperature dependent) 
and expansion rate ( also temperature dependent) of the Universe. When this front factor becomes 
small compared to unity, approach to equilibrium slows down, even if equilibrium is not reached. 
The abundance of the species ip in a comoving volume remains fixed once this factor becomes 
insignificant. This phenomenon is called "freeze out", i.e., the fact that the relative abundance 
does not change after this and continues to evolve like free gas. 

After the species freezes out, at epoch to with corresponding temperature Tjj, the distribution 
function of the species continues to evolve purely due to the effect of the expanding spacetime. 
There are two simple rules of thumb we can prove for its distribution d 3 n/d 3 p in phase space : 

• A relativistic species continues to have the usual Bose-Einstein or Fermi-Dirac distribution 
function (e^ E ± l)" 1 , except that /T 1 = T scales like T(t) = T{t D )R(t D ) / R{t). 

• A species which is non-relativistic, i.e., mass m 3> Tjj the number density simply keeps 
depleting as R~ 3 , just like the particles which are still in equilibrium. But the momenta 
scale as i£ _1 , so energy E — p 2 /2m scales as R~ 2 . This is equivalent to the temperature 
scaling as T(t) = T(t D )R 2 (t D ) / R 2 (t). 

Thus the distribution functions have an energy dependence which is simply obtained from 
their functional forms at the time of decoupling. In the relativistic case in fact remaining self- 
similar, and looks just like that of the particles still in equilibrium, with an important exception. 
If there is a change in the total number of effective degrees of freedom at some temperature, this 
information is not conveyed to the decoupled particles. In the non-relativistic case the scaling of 
the temperature parameter is significantly different. 

9.2 Dark Matter 

There is a variety of evidence to suggest that a large part of the matter content of the Universe is 
neither radiation, nor in the form of baryons. As such it is not capable of participating in processes 
producing electromagnetic radiation and christened Dark Matter. 

The direct evidence for Dark Matter is available at two largely different scales. At the scale 
of individual galaxies and at the scale of clusters of galaxies. At the level of single galaxies it 
is possible to measure speeds of luminous bodies in the spiral arms for those galaxies which are 
visible edge on. The difference in the rcdshifts of the parts rotating away from us and the parts 
rotating towards us is measurable. It turns out that as a function of their distance from the center 
of the galaxy, velocities of rotation in the plane of the galaxy do not slow decrease in accordance 
with the 1/r 2 law expected from Kepler's law. Rather their speeds remain steadily high even 
beyond the visible edge of the galaxy. The plots of the velocity vs. the radial distance from the 
center of the galaxy have come to be called "rotation curves". The departure from Kepler law 
suggests presence of gravitating matter extended to ten times the size of the visible galaxy! 

Secondly at the level of clusters of galaxies, it is possible to measure the relative speeds of 
the galaxies in a cluster, specifically the component of the velocity along the line of sight. By 
viirial theorem the values of these velocities should be set by the total matter content of the 
cluster. Again one finds the velocities more compatible with almost ten times the matter content 
compared to the visible. 
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Another indicator of the extent of the baryonic content is indirect but very sensitive. Big Bang 
Nucleosynthesis predicts ratio of Hydrogen to Helium and the ratios of other light elements to 
Hydrogen determined by one parameter, the baryon to photon ratio, rj = B/s where B is the net 
baryon number (difference of baryon and antibaryon numbers) and the denominator is the photon 
entropy. We shall have occasion to discuss this in greater detail in the section on Baryogenesis. 
The observed ratios of Helium to Hydrogen and other light nuclei to Hydrogen is correctly fitted 
only if 77 ~ 1CT 9 . Knowing the photon temperature very accurately we know the contribution of 
radiation to the general expansion ( it is very insignificant at present epoch). Further knowing 
this accurately we know the baryon abundance rather accurately. Between the two, the latter 
is certainly the dominant contribution to the energy density of the present Universe. However 
the total amount of matter-energy required to explain the current Hubble expansion is almost 30 
times more than the abundance baryons inferred through the BBN data. Again we are led to 
assume the existence of other forms of matter energy that account for the Hubble expansion. It is 
therefore assumed that there is extensive amount of non-relativistic matter present in the Universe 
today, and is called Dark Matter. We do not know at present whether Dark Matter is a single 
species of particles or several different species of particles. We do not know the possible values for 
the masses of such particles, however the study of galaxy formation suggests two classes of Dark 
Matter distinguished by their mass as we see in the next paragraphs. 

The latest data from all sources suggest that the dominant component of the energy driving 
the expansion is actually neither radiation nor matter, but some other form of energy obeying an 
equation of state close to that of relativistic vacuum, p = —p. This is estimated to be contribute 
about 70%. The Dark matter component is estimated to be about 27%, and only about 3% in 
the form of baryonic matter. These conclusions follow from WMAP data on CMB and type la 
supernova data on expansion rate of the Universe more than 7 billion years ago. It is remarkable 
that the approximately 10 times the abundance of Dark Matter relative to baryonic matter as 
inferred directly from galacic and cluster data is verified reasonably accurately by the very indirect 
methods. This is what gives us confidence in the Dark Matter hypothesis. 

When galaxy formation is considered this highly abundant Dark Matter component plays a 
significant role. While no other kind of interaction is permitted between baryonic matter and 
Dark Matter at least at low energies, gravity is always a mediator. It is no surprise therefore that 
the Dark Matter is clustered in approximately the same way as luminous baryonic matter. The 
question whether there are large distributions of Dark Matter separately from baryonic matter 
needs experimentally studied however so far the evidence does not seem to demand such an 
assumption. 

It then follows that the growth of perturbations which led to galaxy formation must have 
proceeded simultaneously for the baryonic matter and the Dark Matter, coupled to each other 
through gravity. The study of this coupled evolution gives rise a distinction of two categories of 
Dark Matter which can be made based on the mass of the corresponding particle. Those particles 
that have become non-relativistic by the time of galaxy formation are called Cold Dark Matter 
( CDM). They are in the form of pressureless dust by this epoch and their chief contribution to 
energy density comes from their rest masses and not their thermal motion, hence Cold. We may 
think of this dividing line as set by the temperature ~ leV when neutral Hydrogen forms. Particles 
which are already non-relativistic at this temperature certainly belong to the category of CDM. 
On the other hand particles that remain a relativistic gas down to leV temperature contribute 
through their thermal energy density and are called Hot Dark Matter ( HDM). A prime candidate 
for this kind of DM is a neutrino, whose masses are constrained to very small values. 

The main difference in the two kinds of DM comes from the nature of the clustering they assist. 
From Jeans formula we see that HDM clustering occurs at large physical scales while CDM can 
cluster at much smaller scales. In fact too much HDM can destroy clumping of baryonic clusters 
at smaller scales. Thus a study of the spectrum of perturbations V(k) gives a clue to the form 
of DM that assisted the formation of galaxies. The current evidence in the light of the WMAP 
data strongly suggests essentially the presence only of CDM, though some proportion of a HDM 
species cannot be ruled out. 

In the following subsections we shall show how we can trace back at least some of the micro- 
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scopic properties of the Dark Matter if we know its abundance today. 



9.2.1 Hot Relics 

For particles that continue to remain relativistic as they are going out of equilibrium, the equa- 
tions from the previous subsection can can be used to show that their abundance at late time is 
determined by the value of their freeze out temperature, i.e., Xf reeze out 

Yoc = Yeq (x freeze out) = 0.278-'^ 



g* s (x) 



If we want to think of this as the Dark Matter candidate, we estimate the energy density it can 
contribute, which is determined to be 



PVo = soYoom = iYoo y-^j keV-(cm) 



From LSS data on distribution of fluctuations, as also the WMAP data it is now concluded that 
the structure formation could not have occurred due to HDM. Hence this is not a very useful 
quantity to verify against observations. Historically, this density value was used to put an upper 
bound on the mass of a neutrino. If the decoupled neutrino is to not be so overabundant that it 
exceeds the current density of the Universe, than its mass must be bounded. 

m < i W »!lM 

9eff 

For j/s the ratio of the g* factors is 0.14, from which one can conclude that m v < 91eV. This is 
known as the Cowsik-McClclland bound. Although the bound is surpassed by both by terrestrial 
experiments and recent astrophysical data, it is an instructive exercise. 

9.2.2 Cold Relics 

For cold relics, we need to determine the quantities Xf, Tf corresponding to the freeze out of the 
species, and its present abundance relative to radiation, Y^. These are determined by solving the 
equation 



It is useful to make an expansion of the cross-section in partial waves, which amounts to an 
expansion in energy, or equivalently in the present setting, an expansion in x = m/T. For a 
massive particle the leading term is 

(a a \v\) = a a (J-^j =cr x- n x > 3 

Thus expressing the cross-section as a function of x, the equation can be solved. The solution to 
this equation gives the left over abundance for a massive particle x at present time. The answer 
typically has the following dependence 



Y x = 0(1) x 



5 / 



( Mp(a A \v\) 



with Xf determined numerically when the Y effectively stops evolving. The present contribution 
to the energy density due to these particles is m x Foo x (s(T ) / p crit ) where s(T ) is the present 
value of entropy density in radiation. 

It is thus possible to relate laboratory properties of the x particle with a cosmological observ- 
able. Given a particle physics model, we can constrain the properties of the potential Dark Matter 
candidate by calculating its contribution to Q dm and then counterchecking the same cross-section 
in collider data. 
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10 Missing Dark Energy 



The important topic of Dark Energy could not be included within the limitations of this course. 
The reader can refer to some of the excellent reviews cited at the end or await the next avatara 
of these notes. 

11 Baryon asymmetry of the Universe 

A very interesting interface of Particle Physics with cosmology is provided by the abundance of 
baryons in the Universe. At first it is a puzzle to note that we only have baryonic matter present 
in the Universe, with no naturally occurring baryons to be seen. 

In principle a cluster of galaxies completely isolated from others could be made totally from 
anti- Hydrogen and anti-elements. However there should be some boundary within which this con- 
fined, since any contact with usual baryonic matter would generate violent gamma ray production 
which would be observed as a part of cosmic rays. But there are no clearly visible empty corri- 
dors separating some galaxies or clusters of galaxies from others, nor is there a significant gamma 
ray background to indicate ongoing baron-anti-baryon annihilation. Thus we assume the present 
Universe to be devoid of priordial anti-baryons. 

Due to charge neutrality the electron number should be exactly equal to the proton number 
of the Universe, and if Lepton number were conserved, we should therefore have a predictable 
abundance of electron type anti-neutrinos. However after the discovery of neutrino oscillations 
the question of total lepton number conservation is open and their number may not be determined 
exactly by the charged lepton number. Thus the total matter vs. anti-matter asymmetry of the 
Universe is a more complicated question. We shall deal only with the baryons where the situation 
is more precisely known. 

The observed asymmetry is quantified by expressing it as a ratio of the photon number density, 
i.e., entropy, 

= ub_ _ n b - n- b 
^ s n 7 

where the upper case subscript B signifies the net baryon number while the lower case subscripts 
b, b signify the actual densities of baryonic and anti-baryonic species separately. Big Bang nu- 
cleosynthesis constrains the value of this ratio very precisely. The abundances of Helium 4 He to 
Hydrogen is sensitively dependent on this ratio, but further, the abundances of light nuclei such 
as Deuterium, 3 He, and 7 Li relative to Hydrogen are also extremely sensitive to this ratio. 

11.1 Genesis of the idea of Baryogenesis 

We believe that the Universe started with a hot Big Bang. If the laws of nature were completely 
symmetric with respect to matter and anti-matter both should be present in exactly same abun- 
dance in thermodynamic equilibrium. Then the asymmetry between the two has to be taken as 
an accidental initial condition. Fortunately we know that the Weak interactions do not respect 
charge conjugation or matter- anti- matter symmetry C, but only the product CP after parity P 
is included. Further, in 1964-65 two crucial discoveries were made. It was shown that certain 
sectors of the theory (K° — K° system) also do not respect CP. In QFT there is a theorem 
that says that it is impossible to write a Lorentz invariant local theory which does not respect 
the combination CPT, now including time reversal T. Thus a CP violating theory presumably 
violates T invariance in same measure. The small mixing of CP eigenstates will also be reflected 
in small asymmetry in reaction rates involving these participants. 

The other crucial discovery was the Cosmic Microwave Background which established the Hub- 
ble expansion back to temperatures as high as 1000K. It is easy to extrapolate this to sufficiently 
early times and higher temperatures when density and temperature would be sufficiently high for 
Particle Physics processes to occur freely. The stage was set for searching for a dynamical expla- 
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nation for the baryon asymmetry (Weinberg 1964) and a specific model was proposed (Sakharov 
1967). 



11.2 Sakharov criteria 

The minimal conditions required for obtaining baryon asymmetry have come to be called Sakharov 
criteria. They can be understood via a specific example. Consider a species X which carries baryon 
number and is decaying into two different possible final products, either two quarks or one anti- 
quark and a lepton. (We use one of the decay modes to determine the baryon number of X and 
violation shows up in the other decay). Such decays are easily possible in Grand Unified models. 
The following should be true for a net B number to remain in the Universe : 



1. Baryon number violation 

X 



qq_ ABi = 2/3 
qi AB 2 = -1/3 



2. Charge conjugation violation 

M(X -> qq) ^ M (X -» qq) 

3. CP violation reflected in difference in rates 

T^X^qq) T^X^qq) _ 

ri= r 1 + r 2 * f 1 + f 2 =ri 

4. Out-of-cquilibrium conditions, which would make reverse reactions become unfavorable 

NetB = ABm + AB 2 (1 - n ) 

(-ABi)fi + (-AB 2 )(l-fi) 

= (ABi - AB 2 ) (n - f x ) 

In the early Universe, the condition for departure from Equilibrium means that the reaction 
rate should become slow enough to be slower than the Hubble expansion rate at that epoch. 
This will happen because reaction cross-sections depend both on density which is falling due to 
expansion, and the energy dependence of the intrinsic cross-section makes it smaller at lower 
temperature. 

T x ^ a x m 2 x /T 
H ~ gl /2 T 2 /M Pl 

1 /2 

Need the rate T x still < H when kT <~ m x . Thus kTjj ~ {a x mpLm x ) 1 . 
Resulting 

~ — x ( Boltzmann evolution) 

s g* 

Thus the result depends purely on the microscopic quantity B (includes Sep) and <?» of the epoch 
when the mechanism operates. 



49 



11.3 Anomalous violation of B + L number 



Quantization of interacting field theories contain many subtleties. Sometimes the process of renor- 
malization does not respect a symmetry present at the classical level. Then quantum mechanically 
the corresponding number operator is not conserved. This situation is called anomaly. This is 
tolerable if it happens for a global charge. If it happens for a gauge charge the model would not 
be viable. It turns out that the Standard Model of Particle Physics does not respect the global 
quantum number B+L, baryon number plus lepton number. The number B — L is also not coupled 
to any gauge charge howeve it remains, miraculously, anomaly free and hence is conserved. 

The anomalous violation is not obtained in perturbation theory. However a handle on the 
anomalous violation rate can be obtained by topological arguments involving global configurations 
of gauge and Higgs fields. A specific configuration whose energy value determines the rate is called 
sphalcron. The energy of a sphaleron is approximately 5 TeV in Standard Model. At temperatures 
well below this, the rate is suppressed exponentially. At a tcmperatrue much higher, the rate is 
order unity. Actually it becomes meaningless to speak of a conserved number. However a number 
generated by any other mechanism will be quickly equilibrated to zero by this non-conservation. 

In the in between regime of temperatures, the rate is estimated as 

r« K (AAV) T 4 e ~E sph (T)/kT 

where k is the determinant of other fluctuations (recall Coleman tunneling formula) and A/Vo 
represents sphaleron zero- mode volume, i.e., the weightage associated with all the possible ways a 
sphaleron can arise. This formula is valid for mw <T< mw/ctw where rmjv is mass of the W 
boson and aw is the fine structure constant <? 2 /47r of the Weak interactions. 

Sphaleron energy depends on the Higgs mass at zero temperature in such a way that too light 
a Higgs (< 90GeV) would result in very rapid violation of B + L around the electroweak phase 
transition. The conclusion is that either the Higgs is heavier ( which is corroborated by the bound 
niH > 117GeV from LEP data), or there is more than one Higgs, or that there was a primordial 
B — L number already present at the electroweak scale. 

11.4 Electroweak baryogenesis 

Could the baryon number arise at the elctroweak scale itself? Sphaleronic processes are already 
becoming unimortant at this scale. Also the properties any new particles needed can be coun- 
terchecked at the LHC or ILC. At the electroweak scale the expansion of the Universe is many 
orders of magnitude (10 12 ) slower than the particle physics processes. Hence direct competition 
with rates is not possible. However, a first order phase transition leads to formation of bubble 
walls. They sweep past any given point in sufficiently short time scale that Particle Physics scales 
compete with this time scale rather than the expansion time scale of the Universe. Such a scenario 
which by-passes the thermal conditions in the Universe is called non-thermal, as against the ex- 
ample studied at the beginning of the section which is called thermal mechanism for baryogenesis. 

If we enhance the SM with additional particles we can actually use the sphaleronic excitations 
to generate B + L asymmetry if the other criteria of Sakharov are satisfied. Typical scenarios rely 
on 

1. Usual C asymmetry of Weak interactions 

2. B + L violation by sphaleronic excitations 

3. CP violation due to complex phases in the vacuum expectation values of one or more scalar 
fields 

4. Out-of-cquilibrium conditions due a first order phase transition 

It turns out that all of these conditions are easily satisfied provided we have more than one Higgs 
scalar and sufficiently large CP phases entering some fermion masses. In specific models favored 
for esthetic reasons however it has not been easy to reconcile all the known constraints from other 
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data with the requirements of clcctroweal baryogenesis. For example, the Minimal Supersymmetric 
Standard Model (MSSM) has the following dangers ( see M. Quiros hep-ph/0101230) 

• Need for first order phase transition implies a light Higgs and a light supcrpartner " stop" of 
the top quark, as also a bound on the ratio of the masses of the two neutral Higgs bosons 
expressed as tan/3, 

110 < m H < llbGeV, tan/3<4, m- tR ~ 105tol65GeV 

• One requires Sep > 0.04 which in turn raises the danger of color breaking vacua. 
11.5 Baryogenesis from Leptogenesis 

A realistic alternative possibility for dynamical explanation for baryon asymmetry is thrown up 
by the discovery of neutrino mass. The very small mass m v ~ 0.01 eV for neutrinos requires 
their Yukawa coupling to the Higgs to be 10" 11 . As we discussed in case of inflation, such small 
dimcnsionlcss numbers seem to hide some unknown dynamics going on. A very elegant explanation 
for the small mass can be given if we assume (i) Majorana type masses for the neutrinos and (ii) 
assume this mass, denoted Mr to be high, Mr <~ 10 14 GeV. It can be shown that 



is a natural relation if such a scale is assumed. Now a scale like 10 14 is also far removed from 
other physics, but is tantalisingly in the range of Grand Unified theories. This mechanism is called 
see-saw mechanism 

This possibility makes leptogenesis naturally possible in the early Universe almost literally 
by the example we studied earlier for the particle X at the beginning of the section. Majorana 
fcrmions do not conserve fermion number. Further, the mixing of the three generations can 
introduce a complex phase in the mass matrix which can lead to CP violation. Finally high 
mass means that the decay rate can actually compete with the expansion scale of the Universe 
which is sufficiently rapid at high temperatures, unlike at electroweak scale. This can result in 
lepton asymmetry of the Universe. This lepton asymmetry converts to baryon asymmetry as 
follows. Recall that at temperatures high compared to the electroweak scale, B + L number is 
meaningless, and will be equilibrated to zero. That is, the anomalous effects ensure A(B + L) = 
and hence will generate approximately AB ~ —AL. The equality is not exact due to interplay of 
several chemical potentials one has to keep track of. 

An important handle on this proposal is provided by the low energy neutrino data. It is possible 
to constrain the extend of CP violation that can be available at high scale from low scale masses 
due to see-saw mechanism. Consider the decay of a heavy neutrino species N into a light lepton 
I and a Higgs particle. There are several such possibilities, and in each case the electric charge in 
the final state is balanced to be zero. Due to lepton number violation characteristic of Majorana 
fcrmions, the same N can also decay into anti-lepton and anti-Higgs. Thus the difference in the 
lepton number of the final products in the two different modes is AL = 2 along the same lines as 
AB = 1 in our example at the beginning of the section. Then the CP-asymmetry parameter in 
the decay of any one of the three heavy neutrinos N t , i — 1,2,3 is defined as 



If we assume a hierarchy of masses Mi < M2 < M3 as is the case of all other fermions, then the 
main contribution to the lepton asymmetry generation comes from the species to decay last, i.e., 
the lightest of the heavy neutrinos Ni. (Why?) The maximum value of CP violation parameter 
d in this case can be shown to be 



m v M R ~ m\ 



w 



T(N t -> £</>) - r(Nj -> l^t) 
T(Ni - U) + T{Ni - ttf) ' 



(3) 




(4) 
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where the mass of the heaviest of the light neutrinos ^3 is bounded by the atmospheric neutrino 
data, which gives the mass-squared difference [S.vr? atm = m\—m\. Thus, m 3 ~ y/ AmJ jm = 0.05eV. 

In the figure [5] we show the solutions of the Boltzmann equations showing the accumulation of 
B — L as temperature T drops for various values of Mi with CP violation chosen to be maximal 
permissible according to above formula and the parameter ra\ = {m D mD)n/Mi chosen 10 _5 eV. 
It turns out that this particular parameter ( numerator is the 11 element of the square of Dirac 
mass matrix for the neutrinos) determines the overall outcome of the thermal B — L number 
production. We see that there is negligible net number B — L at high temperature but it builds 
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Figure 8: The evolution of the B — L asymmetry with temperature, shown here as a function of Z = M\/T, 
with fixed values of Mi as indicated in the legend. The value of the CP violation parameter is maximal 
permissible and the parameter m\ explained in text is chosen 10 eV for all graphs. Figure from N. Sahu 
et al, Nucl. Phy. B752 (2006) 

up as the decay processes are going out of equilibrium. At some point the production levels off. 
Then due to sphalerons, the asymmetry which is initially in the form of light neutrinos also gets 
converted to baryonic form producing net B number. 

From such exercises it can be shown that we need the mass scale Mr of the heavy majorana 
neutrinos to be typically > 10 12 GeV but with some optimism, at least > 10 9 GeV for successful 
thermal leptogenesis. The problem with this conclusion is that firstly a new intermediate scale 
much lower than required for gauge coupling unification is called for. Secondly, as discussed in 
the Introduction, we expect supersymmetry to regulate the QFT involved in Grand Unification 
with several scales of symmetry breaking. But supersymmetry necessarily implies the existence 
of gravitino. Further, it can be shown that if our Universe underwent simple radiation dominated 
expansion from any temperature larger than 10 9 GeV down to Big Bang Nucleosynthesis, sufficient 
number of gravitinos would be generated that would make the Universe matter dominated and 
foul up BBN. Thus it is usual to assume that the "reheat" temperature after inflation is lower 
than 10 9 GeV. But then the thermal leptogenesis discussed here becomes unviable. 

It remains an open possibility that there are non-thermal mechanisms similar to the electroweak 
baryogenesis, but applicable to leptogenesis. 
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12 Appendix 



Here we discuss the "True or False" statements given in section 11.11 Note that some of the 
statements are half baked and warrant a discussion rather than simple yes or no. Some hints. 

1. Curved spacetime takes account of equivalence of gravitational and innertial mass. The 
Relativity principle of space and time could have been Galilean and the formulation would 
be still useful. See ref [1], chapter 12 for Cartan's formulation. 

2. Reparametererization only relabels points. It cannot change physics. Usually the laws are 
written only in the form invariant under rigid rotations. But every law can in principle be 
rewritten to be form invariant under change of parameterization. Thus reparameterization 
invariance cannot be a new physical principle. 

3. Due to Equivalence Principle as adopted by Einstein, all forms of energy are subject to 
and contribute to gravitational field. Energy density therefore must contain contribution of 
gravitational " binding energy" . However we can always choose freely falling frames locally 
so that effect of gravity disappear. In these frames the energy density of gravitational field 
disappears. 

4. Total energy would be an integral of the energy density over a whole spacelike surface. 
This answer would remain unchanged under coordinate local transformations especially if 
we restrict ourselves to rigid transformations at infinity ( sitting where we measure up the 
energy). But GR throws up the possibility of compact spacelike hypersurfaces. In this case 
asymptotic region is not available. 

5. If this genuinely means spacetime measurements are meaningless at that point then it is 
unphysical. But it can be an artifact of coordinate system, as for instance the origin in a 
spherical or cylindrical coordinates. 

6. Divergence of metric coefficients is often avoided using different coordinate systems. 

7. Curvature tensor is a physical quantity. Divergence of its components will also often imply 
divergence of some components of energy-momentum tensor. Such points would be unphys- 
ical. However note that much electrostatics is done assuming point charges. These have 
infinite energy density at the location of the point. When such points are isolated we hope 
some other physics takes over as the singular point is approached. 

8. The expansion of the Universe is neither relativistic, nor a strong gravity phenomenon at 
least ever since BBN. It admits Newtonian description. If the spacelike hypersurfaces were 
compact that would be easier to explain as a dynamical fact in GR. In Newtonian physics 
we would simply accept is as fact, just as we are willing to accept infinite space as fact. 
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